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ABSTRACT 

We calculate accurate solar models and report the detailed time dependences of 
important solar quantities. We use helioseismology to constrain the luminosity evolution 
of the sun and report the discovery of semi-convection in evolved solar models that 
include diffusion. In addition, we compare the computed sound speeds with the results 
of p-mode observations by BiSON, GOLF, GONG, LOWL, and MDI instruments. We 
contrast the neutrino predictions from a set of eight standard-like solar models and 
four deviant (or deficient) solar models with the results of solar neutrino experiments. 
For solar neutrino and for helioseismological applications, we present present-epoch 
numerical tabulations of characteristics of the standard solar model as a function of 
solar radius, including the principal physical and composition variables, sound speeds, 
neutrino fluxes, and functions needed for calculating solar neutrino oscillations. 
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1. Introduction 

Why are new calculations of standard solar models of interest? After all, solar models have 
been used to calculate neutrino fluxes since 1962 (Bahcall, Fowler, Iben, & Sears 1963) and solar 
atmospheres have been used to calculate p-mode oscillation frequencies since 1970 (Ulrich 1970; 
Leibacher & Stein 1971). Over the past four decades, the accuracy with which solar models are 
calculated has been steadily refined as the result of increased observational and experimental infor- 
mation about the input parameters (such as nuclear reaction rates and the surface of abundances of 
different elements), more accurate calculations of constituent quantities (such as radiative opacity 
and equation of state), the inclusion of new physical effects (such as element diffusion), and the 
development of faster computers and more precise stellar evolution codes. 

Solar models nevertheless remain at the frontiers of two different scientific disciplines, solar neu- 
trino studies and helioseismology. In an era in which many major laboratory studies are underway 
to study neutrino oscillations with the aid of very long baselines, ~ 10 3 km, between accelerator and 
detector, solar neutrinos have a natural advantage, with a baseline of 10 8 km (Pontecorvo 1968). 
In addition, solar neutrinos provide unique opportunities for studying the effects of matter upon 
neutrino propagation, the so-called MSW effect (Wolfenstein 1978; Mikheyev & Smirnov 1985), 
since on their way to terrestrial detectors they pass through large amounts of matter in the sun 
and, at night, also in the earth . 

The connection with ongoing solar neutrino research imposes special requirements on authors 
carrying out the most detailed solar modeling. Precision comparisons between neutrino measure- 
ments and solar predictions are used by many physicists to refine the determination of neutrino 
parameters and to test different models of neutrino propagation. Since the neutrino experiments 
and the associated analysis of solar neutrino data are refined at frequent intervals, it is appropriate 
to reevaluate and refine the solar model predictions as improvements are made in the model input 
parameters, calculational techniques, and descriptions of the microscopic and macroscopic physics. 

In this paper, we provide new information about the total solar neutrino fluxes and the pre- 
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dicted neutrino event rates for a set of standard and non-standard solar models. Using the best- 
available standard solar model, we also present the calculated radial dependence of the production 
rate for each of the important solar neutrino fluxes. We publish for the first time the results of a 
precision calculation with the standard solar model of the electron density throughout the sun, from 
the innermost regions of the solar core to the solar atmosphere. We also present for the first time 
a detailed calculation of the radial profile of the number density of scatterers of sterile neutrinos. 
These quantities are important for precision studies of neutrino oscillations using solar neutrinos. 

We also provide detailed predictions for the time evolution of some of the important solar 
characteristics such as the depth and mass of the solar convective zone, the radius and the luminosity 
of the sun, the central temperature, density, pressure, and hydrogen mass fraction, as well as the 
temperature, density, pressure, and radiative opacity at the base of the convective zone. As far 
as we know, these are the first detailed results submitted for publication on the time evolution of 
many of these quantities. Some of the calculated time dependences may be subject to observational 
tests. 

At the present writing, the sun remains the only main-sequence star for which p-mode os- 
cillations have been robustly detected. Thus only for the sun can one measure precisely tens of 
thousands of the eigenfrequencies for stellar pressure oscillations. The comparison between the 
sound speeds and pressures derived from the observed p-mode frequencies and those calculated 
with standard solar models has provided a host of accurate measurements of the interior of the 
nearest star. The solar quantities determined by helioseismology include the sound velocity and 
density as a function of solar radius, the depth of the convective zone, the interior rotation rate, 
and the surface helium abundance. The excellent agreement between the helioseismological ob- 
servations and the solar model calculations has shown that the large discrepancies between solar 
neutrino measurements and solar model calculations cannot be due to errors in the solar models 
(cf. Figure 11). 

In this paper, we present a refined comparison between our best standard solar model and 
measurements of the solar sound speeds obtained using oscillation data from a number of different 
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sources. We use a combination of data from the LOWL instrument and the BiSON network, two 
sets of data from the GOLF instrument, as well as data from the GONG network and the MDI 
instrument. 

We describe in the following paragraph the organization of this paper. But, since this paper 
contains a lot of information on disparate topics, we recommend that the reader first turn to the 
final section, § 8, and peruse our summary and discussion of the main new results. Section 8 may 
help the reader decide which sections of the paper the reader wants to read (or to skip). The 
different sections are written so that they can be read more or less independently. 

Section 2 defines the slightly-improved Standard solar model and describes the numerical 
tables that present details of the contemporary characteristics of the Standard model. Section 3 
describes the main sequence time dependences of some of the most important characteristics of the 
Standard model. We present in this section the time evolution of the solar radius and luminosity, 
the properties of the convective zone, and the physical characteristics of the center of the sun. We 
discuss solar properties from the initial age main sequence to an age of 8 billion years. Section 4 
defines and compares the physical characteristics of seven variant standard models and four deviant 
(deficient) solar models, which together with the Standard model make up a set of twelve models 
whose neutrino fluxes we evaluate in this paper. We previously used a subset of nine of these 
models to test the robustness of helioseismological inversions (Basu, Pinsonneault, & Bahcall 2000, 
hereafter BPB2000). We discuss in this same section two new standard-like solar models with heavy 
element to hydrogen ratios that differ slightly from our previously adopted value of Z/X. Section 5 
discusses solar neutrino physics. We present the predicted present-day neutrino fluxes for the 
Standard model and for all of the variant and deviant solar models, as well as the electron number 
density versus solar radius. We contrast the predicted neutrino event rates with the results of the 
chlorine, Kamiokande, GALLEX, SAGE, Super-Kamiokande, and GNO solar neutrino experiments. 
We also give in this section the calculated time evolution of the most important solar neutrino 
fluxes. Section 6 compares calculated and observed sound speeds. We present the results both 
on a panoramic scale suitable for discussing the implications for solar neutrino physics and on a 



- 5 - 



zoomed-in scale appropriate for detailed investigations of solar physics. In § 7, we compare, for all 
12 of the solar models discussed in this paper, the calculated values of the surface helium abundance 
and the depth of the convective zone with the measured values for these quantities. We summarize 
and discuss our main results in § 8. 

The interested reader may wish to consult the following works that summarize the solar neu- 
trino aspects of solar models (Bahcall 1989; Bahcall & Pinsonneault 1992, 1995; Berezinsky, Fioren- 
tini, & Lissia 1996; Castellani et al. 1997; Richards et al. 1996; Turck-Chieze et al. 1993; Bahcall, 
Basu, & Pinsonneault 1998) and the helioseismologic aspects of solar models (Bahcall & Ulrich 
1988; Bahcall & Pinsonneault 1995; Christensen-Dalsgaard et al. 1996;; Guenther & Demarque 
1997; Guzik 1998; Turck-Chieze et al. 1998; Brun, Turck-Chieze, & Zahn 1999; Fiorentini & Ricci 
2000). 

2. Standard solar model 

We define in § 2.1 the Standard solar model and present in § 2.2 some of the important 
contemporary characteristics of the Standard solar model, including detailed tables of the physical 
variables as a function of the solar radius. 

By 'the Standard solar model', we mean the solar model which is constructed with the best- 
available physics and input data. All of the solar models we consider, standard or 'deviant' models, 
(see below) are required to fit the observed luminosity and radius of the sun at the present epoch, as 
well as the observed heavy element to hydrogen ratio at the surface of the sun. No helioseismological 
constraints are used in defining the Standard solar model. 

Naturally, Standard models improve with time, as the input data are made more accurate, 
the calculational techniques become faster and more precise, and the physical description is more 
detailed. Thus it is necessary at each presentation of a new Standard model to describe what has 
changed from the previous best model and to provide references to the earlier work. The reader 
can see Bahcall (1989) for a general reference to the early work on precise solar models that were 
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motivated by the need to predict accurate solar neutrino fluxes and to make comparisons with 
helioseismological data. 

2.1. Definition of the Standard model 

Our standard solar model 1 is constructed with the OPAL equation of state (Rogers, Swenson, 
& Iglesias 1996) and OPAL opacities (Iglesias k, Rogers 1996), which are supplemented by the low 
temperature opacities of Alexander & Ferguson (1994). The model was calculated using the usual 
mixing length formalism to determine the convective flux. 

The principal change in the input data is the use of the Grevesse & Sauval (1998) improved stan- 
dard solar composition in the OPAL opacities (see http://www-phys.llnl.gov/Research/OPAL/index.htm) 
and in the calculation of the nuclear reaction rates. The refinements in this composition redeter- 
mination come from many different sources, including the use of more accurate atomic transition 
probabilities in interpreting solar spectra. The OPAL equation of state and the Alexander and Fer- 
guson opacities are not yet available with the composition recommended by Grevesse and Sauval 
1998. 

We have used a present-epoch solar luminosity of 1369 Wm 2 [47r x (A.U.) 2 ] for all the models 
described in detail in this paper. Only for the neutrino fluxes presented in § 5.1.1 and § 5.1.3 have we 
used the more recent best-estimate value of solar luminosity, 1366.2 Wm 2 [4-7r x (A.U.) 2 ] = 3.842 x 
10 33 ergs -1 (see Frohlich & Lean 1998; Crommelynck, Fichot, Domingo, & Lee 1996). The differ- 
ence between these two values for the luminosity is 0.2%. For the calculations of uncertainties in 
neutrino flux predictions, we assume a la uncertainty of 0.4%. The uncertainty in the predicted 
solar neutrino fluxes due to the luminosity is an order of magnitude smaller than the total uncer- 
tainty in predicting the neutrino fluxes. For all the other quantities we calculate, the uncertainty 

lr To simplify the language of the discussion, we will often describe characteristics of the Standard model as if we 
knew they were characteristics of the sun. We will sometimes abbreviate the reference to this Standard model as 
BP2000 or Bahcall-Pinsonneault 2000. 
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in the solar luminosity has an even smaller effect. 

The nuclear reaction rates were evaluated with the subroutine exportenergy.f (cf. Bahcall &; 
Pinsonneault 1992), using the reaction data in Adelberger et al. (1998) and with electron and 
ion weak screening as indicated by recent calculations of Gruzinov & Bahcall (1998); see also 
Salpeter (1954) 2 . The model incorporates helium and heavy element diffusion using the exportable 
diffusion subroutine of Thoul (cf. Thoul, Bahcall & Loeb, 1994; Bahcall & Pinsonneault 1995) 3 . 
An independent and detailed treatment of diffusion by Turcotte et al. (1998b) yields results for the 
impact of diffusion on the computed solar quantities that are very similar to those obtained here. 
We have used the most recent and detailed calculation (Marcucci et al. 2000a) for the So-factor 
for the 3 He(p,e + + ^ e ) 4 He reaction: So (hep) = 10.1 x 10~ 20 keV b, which is a factor of 4.4 times 
larger than the previous best-estimate [see § 5.1.2, Bahcall & Krastev (1998), and Marcucci et al. 
(2000b) for a discussion of the large uncertainties in calculating So (hep)]. For values of So (hep) in 
the range of current estimates, the assumed rate of the hep reaction only affects in a noticeable way 
the calculated flux of hep neutrinos and does not affect the calculated fluxes of other neutrinos, the 
helioseismological characteristics, or other physical properties of the sun. 

For the standard model, the evolutionary calculations were started at the main-sequence stage. 
The model has a radius of 695.98 Mm. We do not study the pre-main sequence evolution in our 
standard model. This epoch has been discussed in the classical work of Iben (1965) and the effects 
of pre-main sequence have been shown by Bahcall & Glasner (1994) to be unimportant for solar 
neutrino calculations (see also the thorough discussion by Morel, Provost, & Berthomieu 2000). We 

Other approximations to screening are sometimes used. The numerical procedures of Dzitko et al. (1995) and 
Mitler (1977) predict reaction rates that are too slow for heavy ions because they assumed that the electron charge 
density near a screened nucleus is the unperturbed value, en e (oo). This assumption seriously underestimates the 
charge density near heavy ions. For example, it is known that a screened beryllium nucleus under solar interior 
conditions has charge density near the nucleus ss — 3.85en e (oo) (Gruzinov & Bahcall 1997; Brown & Sawyer 1997; all 
quantum mechanical calculations give similar results, see Bahcall 1962, and Iben, Kalata, & Schwartz 1967). 

3 Both the nuclear energy generation subroutine, exportenergy.f, and the diffusion subroutine, diffusion. f, arc 
available at the Web site www.sns.ias.edu/~jnb, menu item: neutrino software and data. 
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do consider one pre-main sequence model, which differs very little from the corresponding model 
started at the zero-age main sequence. 

The ratio of heavy elements to hydrogen (Z/X) at the surface of the model is 0.0230, which was 
chosen to be consistent with the value obtained by Grevesse & Sauval (1998). A Krishna-Swamy T- 
t relationship for the atmosphere was used. We adopt a solar luminosity L Q = 3.844 x 10 33 erg s" 1 
and a solar age of 4.57 x 10 9 yr (see Bahcall & Pinsonneault 1995). 

In the course of our recent analysis of systematic uncertainties in helioseismological inversions 
(Basu, Pinsonneault, & Bahcall 2000), we uncovered an error in the code we wrote for the interpo- 
lation of the OPAL 95 opacities. The edges of opacity tables were flagged by zero opacity values; 
unfortunately, there were some interpolation problems associated with the occurrence of zero values 
inside the table. This problem occurred because the logarithm of the opacity, which is what we 
were tabulating, can actually be zero. The interpolation error caused small changes in the opacity 
which produced errors in the sound speed of order 0.1% for solar radii in the range of O.3i?0 to 
0.7 Rq and errors of order 3% in the 8 B neutrino flux (with smaller errors for other neutrino fluxes). 
In the following sections, we will point out more specifically the changes in the neutrino predictions 
and in the sound velocities that are produced by correcting this interpolation error. These changes 
are particularly apparent in the comparison of Figure 13 and Figure 14, which are discussed in the 
helioseismology section, § 6. 



2.2. Some contemporary characteristics of the standard solar model 

The details of the structure of a standard solar model are of importance for both helioseismol- 
ogy and solar neutrino calculations. In the first paper in this series in which we published the details 
of a model of the solar interior structure (Bahcall et al. 1982), we presented a table with only 27 
radial shells (rows of numbers) and 10 variables for each shell (mass, radius, temperature, density, 
hydrogen fraction, helium fraction, luminosity generated, and the source density of the p-p, 7 Be, 
and 8 B neutrinos). Over the years, much greater precision was required by the increasing sophisti- 
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cation of the helioseismological measurements and the solar neutrino calculations. Fortunately, the 
computing capacity more than made up for the necessary increase in model details. 

We have created from the output of the present calculations two exportable computer files that 
contain in easily readable form the details of our best Standard solar model (BP2000). Physical 
variables are given at 875 separate radial shells, which is the same number of shells used to calculate 
the solar interior model. In addition to the variables cited above, this file contains the pressure, 
electron number density, the mass fractions of 3 He, 7 Be, 12 C, 14 N, and 16 0, as well as the source 
densities of all eight of the most important solar neutrino fluxes. These files are accessible at 
http://www.sns.ias.edu/~jnb. Previous standard solar models in this series (published in 1982, 
1988, 1992, 1995, as well as 1998) are available at the same URL and can be used to test the 
robustness of conclusions that depend upon models of the solar interior. 



3. Time dependences 

In this section, we present and discuss the time dependences of some of the principal character- 
istics of the standard solar model. We begin in § 3.1 by describing the separate temporal evolution 
of the solar radius and the solar luminosity and then discuss the simple relation we have found 
between the values of these two quantities as a function of time. We also demonstrate that the 
solar luminosity is a robust function of time. In § 3.2, we present the time-dependent fractions of 
the solar luminosity that are produced by different nuclear fusion reactions. We concentrate in § 3.3 
on the convective zone, presenting the calculated time dependences of the depth and the mass of 
the convective zone. We also report the time dependence of the temperature, density, pressure, and 
opacity at the base of the convective zone. In § 3.4, we discuss the time dependence of quantities 
at the center of the sun, the central temperature, density, pressure, and hydrogen mass fraction. 
In § 3.5, we calculate and discuss the large and small separations of the p-mode frequencies as a 
function of solar age. Since we do not discuss pre-main sequence evolution, our calculations are not 
precise for times less than 0.1 x 10 9 yr (see Morel, Provost, & Berthomieu 2000; Weiss & Schlattl 
1998). 
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Iben (1967,1974) and Demarque & Guenther (1991) summarize in comprehensive presentations 
the evolution of solar parameters in models that were calculated prior to the inclusion of element 
diffusion in solar evolutionary codes. These discussions did not encounter the problem of semi- 
convection discussed here in § 3.3 because this phenomenon is caused by effects of diffusion near 
the base of the convective zone. 

The solar radius and luminosity (or equivalently, the solar effective temperature and luminos- 
ity) constitute precise constraints on the possible geological histories of the earth. We quantify 
these constraints in the following subsection and specify upper limits to the allowed discrepancies 
from the standard solar model profile of solar luminosity versus age. 
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Fig. 1. — The calculated radius, Ro(t), as a function of age for the standard solar model, Bahcall- 
Pinsonneault (2000). The solar age is measured in units of 10 9 yr. The present age of the sun, 
4.57 x 10 9 years, is indicated by an arrow in Figure 1. The radius increases from 0.87-Rq at the 
zero age main sequence to 1.0-R© at the present epoch and 1.18-R© at a solar age of 8 billion years. 
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Table 1: The solar radius, R®(t), as a function of solar age for the Standard solar model. The 
present age of the sun is 4.57 x 10 9 years (see Bahcall k, Pinsonneault 1995) 
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3.1. Radius and luminosity versus age 



Figure 1 and Table 1 present the calculated radius of the sun versus solar age for the standard 
solar model. The values given are normalized to the present-day solar radius, Rq. Over the lifetime 
of the sun, the model radius has increased monotonically from an initial value of O.869i?0 to the 
current value 1.0i? Q , a 15% rise. At a solar age of 8 billion years, the solar radius will increase 
to 1.17i? Q . We shall see in the following discussion that some of the important evolutionary 
characteristics of the sun can be expressed simply in terms of the ratio R & (t)/R Q (today). 

Figure 2 and Table 2 show the calculated solar luminosity as a function of solar age, normalized 
to the present-day solar luminosity, L Q . For the standard model, the total solar luminosity has 
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Fig. 2. — The normalized solar luminosity, L@(t)/L@ (today) versus solar age for the Standard solar 
model (solid curve) and for three 'deficient' solar models: the No Diffusion model (dotted curve), 
the S34 = model (short dashes), and the Mixed model (long dashes). The luminosity evolution 
of the sun is essentially the same in all solar models we have investigated, including deficient solar 
models. The rms deviation of the deviant models from the standard solar model luminosity is only 
1% over the history of the sun from the zero- age main sequence to the current epoch (see text for 
more details). The product L@(i) x R & (t)~ 2 ' 5 varies by ±4% over the entire period from the zero 
age main sequence to a solar age of 8 billion years, while the solar luminosity itself varies by slightly 
more than a factor of two during this period. In the period between 4 billion years to 8 billion 
years, the relation L®(t) oc -R©(i) 2 is satisfied to ±1/2%. The solar luminosity has increased by 
48% from the zero main sequence to the present epoch. The present age of the sun is indicated by 
an arrow at 4.59 x 10 9 years. 

risen monotonically from a zero- age value of 0. 6771*0 . 

The time evolution of the solar luminosity is robust. We also show in Figure 2 the solar 
luminosity as a function of time for the three most deficient solar models that are described in 
the following section, § 4. The rms difference between the standard luminosity and the luminosity 
of the deviant models is 1.6% for the mixed model (1.2% ignoring the first Gyr), 0.7% for the no 
diffusion model (0.5% ignoring the first Gyr), and 0.9% for the S34 = model (0.8% ignoring the 
first Gyr). The largest deviations occur for the zero-age main sequence models and are 2.5% for 
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Table 2: Solar luminosity as a function of solar age for the Standard solar model. The table gives 
the computed values of the solar luminosity in units of the present-day solar luminosity. The model 
was iterated to give the observed luminosity at the present epoch, 4.57 x 10 9 yr. 
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the mixed model, 1.9% for the no diffusion model, and 1.7% for the = model. All of the solar 
models show essentially the same shape for the luminosity evolution as a function of age. 

Figure 2 shows that the product 

L Q (t) x R Q (t)~ 2 - 5 = constant (1) 

to an accuracy of about ±4% over the entire period from the zero age main-sequence to an age 
of 8 billion years. The solar luminosity L©(i) varies from O.68L to L at the present epoch 
and will reach 1.36L after 8 billion years on the main sequence. The corresponding values of 
L (i) x [i? (today)/i? (t)] 2 - 5 are 0.96, 1.00, and 0.93. Thus the luminosity of the sun varies by 
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slightly more than a factor of two while L©(i) x [i?Q(today)/i? Q (t)] 2 ' 5 varies by only a few percent. 

Over the first 5 billion years, Eq. (1) is satisfied to an accuracy of ±2%. From 4 billion years 
to 8 billion years, the relation is somewhat less steep, L Q (t) oc R&(t) 2 to an accuracy of ±1/2%. 
This transition from L (£) oc R & (t) 2 - 5 to L & (t) oc -R©(t) 2 can be seen clearly in Figure 2. 

Table 3: The calculated effective temperature as a function of solar age for the Standard model. 
The solar age is measured in units of 10 9 yr and the solar effective temperature is measured in units 
of the present-epoch effective temperature 
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Table 3 gives the calculated effective temperature of the Standard model, T c g, as a function 
of solar age. The effective temperature varies by only ±1.5% over the entire period from the zero 
age main sequence to an age of 8 billion years and only by ±0.7% from 2 billion years to 8 billion 
years. 

The relation between L Q (t) and i?©(i) discussed above can be restated using the calculated 
time-dependence of the solar effective temperature, which is summarized in Table 3. In the period 
between 4 Gyr and 8 Gyr, the effective temperature is essentially constant (to an accuracy of 
±0.25%). This constancy of the effective temperature results in the scaling relation L Q (t) oc i?©(i) 2 
that is valid in this period. The effective temperature changes by an order of magnitude larger 
fraction during the evolution up to the present-age sun, which results in a dependence closer to 
L & oc Rq 5 (cf. also Figure 2 of Demarque & Guenther 1991). 
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Fig. 3. — The fraction in the Standard model of the solar luminosity produced by different nuclear 
fusion reactions versus solar age. The luminosity generated by the p — p nuclear fusion branch 
that is terminated by the 3 He- 3 He reaction is marked by a dashed curve in the figure and the 
luminosity produced by the p — p branches that proceed through the 3 He- 4 He reaction is denoted 
by a dot-dashed curve. The luminosity generation by the CNO cycle is indicated by a solid line. 
The unit of luminosity is the present-day total solar luminosity. At the present epoch, the p — p 
+ 3 He- 3 He reactions produce 87.8% of the solar luminosity and the branches terminating through 
the 3 He- 4 He reaction generate 10.7% of the solar luminosity. The CNO cycle produces 1.5% of the 
present-epoch luminosity. 



3.2. Energy fractions 



Figure 3 shows, for the Standard model, the energy generated by different nuclear fusion 
reactions as a function of solar age. The present-day total solar luminosity, Lq (today), is the unit 
of luminosity in Figure 3. 

The branch of the p — p chain that is denoted in Figure 3 by p — p + 3 He- 3 He (the dashed 
curve) proceeds primarily through the reactions p(p, e + i/ e ) 2 H(j>, 7) 3 He( 3 He, 2p) 4 He. For simplic- 
ity, we include all p — p reactions in this sum but do not show explicitly the pep reactions in 
the above scheme. The small energy contribution due to pep reactions is included in the calcu- 
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lations that led to Figure 3. The 3 He- 4 He branch (the dot-dashed curve) includes the nuclear 
reactions that produce both the 7 Be and the 8 B neutrinos: 3 He( 4 He, 7) 7 Be(e~, u e ) 7 Li(p, 4 He) 4 He 
and 3 He( 4 He, >y) 7 Be(p, 7) 8 B( 4 He+e++^ e ) 4 He. The CNO reactions are denoted by CNO in Figure 3. 

The branch that terminates via the 3 He- 3 He reaction dominates the solar energy generation 
throughout the main sequence lifetime shown in Figure 3, although the CNO reactions overtake 
the 3 He- 4 He branches at an age of about eight billion years. At an age of one billion years, 96.7% 
of the solar luminosity is generated through p — p reactions plus the 3 He- 3 He termination, 2.8% 
through the 3 He- 4 He termination, and only 0.5% through the CNO cycle. The situation gradually 
changes as the sun heats up. At the present epoch, the 3 He- 3 He termination produces 87.8% of the 
solar luminosity and the branches terminating through the 3 He- 4 He reaction generate 10.7% of the 
solar luminosity. The CNO cycle produces 1.5% of the present-epoch luminosity. By the time the 
Standard solar model reaches an age of eight billion years, the percentages are 57.6%, 20.4%, and 
22.0%, respectively, for 3 He- 3 He, 3 He- 4 He, and CNO reactions. The energy loss due to gravitational 
energy expansion ranges from —0.03% to —0.04% to —0.07% as the sun's age increases from one 
billion years to the present epoch to eight billion years. 



3.3. Convective zone 

Figure 4 shows the depth of the convective zone,i?(CZ, t), as a function of age for the standard 
solar model. Correspondingly, Figure 5 shows the mass of the convective zone, M(CZ,t), as a 
function of age. In both cases, the temporal dependence from the zero age main sequence to the 
current epoch is describable by a simple function of Ro(t)/R@ (today). 

We limit the time period covered in Figure 4 and Figure 5 to be less than 6.5 billion years, 
since between 6.5 and 7.0 billion years of age semi-convection begins to influence the computed 
R(CZ,t) and M(CZ,i). The evolution of the depth of the solar convective zone was previously 
studied by Demarque & Guenther (1991) in an investigation that did not include element diffusion. 

The occurrence of semi-convection in evolved solar models is discussed here for the first time, 
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Fig. 4. — The depth of the convective zone, i?(CZ, t), as a function of age for the standard solar 
model. The depth of the convective zone is approximately proportional to the contemporary solar 
radius. The solar age is measured in units of 10 9 yr. 

to the best of our knowledge. The onset of semi-convection is triggered by the effects of element 
diffusion, which was first included in detailed solar models in the early 1990's (see, e.g., Bahcall & 
Loeb 1990; Proffitt & Michaud 1991; Bahcall & Pinsonneault 1992). 

The process works essentially as follows. The time scale for diffusion decreases as the surface 
convective zone becomes shallower; the metal abundance increases steadily below the surface con- 
vective zone. At some point in the evolutionary history, near an age of 6.5 billion years, the opacity 
from the enhanced metal abundance below the surface convective zone becomes large enough to 
make the metal-rich radiative layers below the surface convective zone convectively unstable. How- 
ever, the mixing of the metals into the envelope causes a local drop in the metal abundance Z and 
and the opacity k at the base of the convective zone. This result in turn causes the convective zone 
to recede until the metal abundance builds up again. Richer, Michaud, & Turcotte (2000) discuss 
a related phenomenon in A and F stars. 
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We noticed the existence of semi-convection in our standard solar model only because we made 
precise plots of the depth and mass of the convective zone as a function of time (see Figure 4 and 
Figure 5). The effects of semi-convection were not noticeable in plots of external quantities such 
as the solar luminosity or effective temperature (see § 3.1 ). In our current code, if a region is 
convectively unstable according to the Schwarzschild criterion, then we instantly mix this material 
with the other material in the convective zone. A more accurate treatment, which allows for the 
possibility of a composition gradient being established at the base of the convective zone, is required 
to calculate reliably the influence of semi-convection once it begins to be perceptible in the numerical 
solutions [see, for example, Merryfield (1995), Canuto (2000) and references quoted therein]. This 
is why we have terminated the plots ofi?(CZ, t) and M(CZ, t) at a solar age of 6.5 x 10 9 yr. 
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Fig. 5. — The mass included within the convective zone, M(CZ,t), as a function of age for the 
standard solar model. The mass of the convective zone is measured in units of Mq (today) and 
the solar age is measured in units of 10 9 yr. The mass of the convective zone is approximately 
proportional to Ro(t)~~ 2 . 



From the zero age main sequence to the present solar age, the depth of the convective zone is 
approximately proportional to the contemporary solar radius, i?©(t), i.e., 
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0.714 < R(CZ,t)/R Q (t) < 0.722. 



(2) 



The zero-age main sequence value is R(CZ, t = 0) = 0.627i?o (today) and the present-day value is 
R(CZ, today) = 0.714R Q (today), which corresponds to an approximately 14% decrease in the depth 
of the convective zone over the age of the sun. However, if we normalize the depth of the convective 
zone to the contemporary solar radius, the change is very small, R(CZ,t = 0) = 0.722R & (t = 0) 
and R(CZ, today) = 0.714i? Q (today). After 6.5 billion years, R(CZ,t = 6.5) = O.764i? (today) 
and, in terms of the contemporary solar radius, R(CZ,t = 6.5) = O.7O6i? (t = 6.5). The relative 
evolution between the convective zone depth and the total solar radius, d[R(CZ, t)/R®(t)]/dt = 
—0.002 per 10 9 yr, is very slow. 

Over the period from the initial age main sequence to the present epoch, the mass of the solar 
convective zone decreases with time as R & (t)~ 2 , 



to an accuracy of about ±1%. The zero-age main sequence value of the mass included within the 
convective zone is 0.0329M Q , which decreases to 0.02415M Q at the present epoch. However, this 
proportionality is no longer valid for larger ages, as can be seen in Figure 5. The mass of the 
convective zone is M(CZ, 6.5) = 0.0233M© at 6.5 billion years. At about 6.8 billion years, the 
previous monotonic decrease of the mass of the convective zone is reversed and M(CZ,t) begins 
to increase with time. This behavior is not shown in Figure 5 since the calculated increase in the 
mass of the convective zone occurs in the same time frame as the onset of semi-convection. 

Figure 6 shows some properties at the base of the convective zone as a function of solar age. 
The figure displays the time dependence of the temperature, T(CZ,t) (solid curve), the density, 
p(CZ,i) (dotted), the pressure, P(CZ,t) (short dashes), and the radiative opacity, re(CZ,i) (long 
dashes). For convenience in plotting, each of the physical variables has been divided by its value at 
the current epoch. In cgs units, the Standard model parameters have the following values at an age 
of 4.57 x 10 9 yr: T(CZ, today) = 2.18 x 10 6 K, p(CZ, today) = 0.19, P(CZ, today) = 5.58 x 10 13 , 



M(CZ,t) x R G (t) 2 = constant 



(3) 
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Fig. 6. — Some properties at the base of the convective zone as a function of age for the Standard 
solar model. The properties shown are the temperature, T(CZ,t) (solid curve), the density, p(CZ, t) 
(dotted), the pressure, P(CZ,t) (short dashes), and the radiative opacity, n(CZ,t) (long dashes). 
All of the quantities are normalized by dividing by their values at the present epoch. After 6.5 
billion yr, semi-convection begins to be important. 

and k(CZ, today) = 20.5. 

The temperature at the base of the convective zone decreases by 14% from the zero-age main- 
sequence to the current solar age, i. e., to a good approximation 



T(CZ,i) oc RiC^ty 1 . (4) 

The opacity at the base of the convective zone decreases by 6% over the same period. The density 
and the pressure at the base of the convective zone decrease by much larger quantities, by factors of 
2.1 and 2.3, respectively. Equation 4 is valid to an accuracy of about ±1% throughout the 6.5 x 10 9 
yr shown in Figure 6. 

The base of the convective zone is defined by the Schwarzschild criterion. Because the adiabatic 
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gradient is nearly constant in time, the equations of stellar structure imply that the quantity 
kPL/MT 4 will be approximately constant at the base of the convective zone. From Figure 6, we 
see that both the opacity and the temperature decrease slowly at the base of the convective zone. 
Solar models therefore compensate for the increase of the luminosity by the decrease of the pressure 
at the boundary between radiative and convective equilibrium. 



3.4. Central values of Temperature, Density, and Pressure 

Figure 7 shows the time dependence of the central values for the temperature, density, and 
pressure of the Standard solar model. The results are normalized to the computed values for the 
present epoch. 
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Fig. 7. — The temporal evolution of the central temperature, density, pressure, and hydrogen 
mass fraction. The figure shows the computed values for the Standard solar model of the central 
temperature (solid line), pressure (dot-dash line), density (dash line), and hydrogen mass fraction 
(dotted line). 



Over the 8 billion years shown in Figure 7, the central temperature increases by about 39%. 
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At the zero-age main sequence, the central temperature is 13.5 x 10 9 K and is 15.7 x 10 9 K at the 
current epoch; the central temperature of the model reaches 18.8 x 10 9 K at a solar age of 8 billion 
years. This increase is very similar to that of the solar radius. In fact, the central value of the 
temperature is, to a reasonable approximation, proportional to the solar radius, i?©(i). The ratio 



to an accuracy of ±1.5% (±3%) from the zero ago-main sequence to the current epoch (to a solar 
age of 6.5 billion years). 

The sun derives its luminosity by burning hydrogen into helium. The hydrogen mass fraction 
at the center of the sun, X c , decreases from 0.708 at the zero-age main sequence to 0.340 at the 
current epoch, a decrease by more than a factor of two. At an of 6.5 billion years, the mass fraction 
has decreased to 0.145, a decrease by a factor of almost five from the zero-age value. The remaining 
hydrogen is burned very quickly, with X c reaching 0.008 at a solar age of 8 billion years. 



The p-mode oscillation frequencies change as the structure of the sun evolves with age. Since 
only low degree modes (small spherical harmonic, I) are observable for stars other than the Sun, 
we shall discuss the changes in just these modes. 

For the high-order (large radial node, n), low-degree modes that are expected to be observed 
in stars, the frequencies satisfy the following relation to a good approximation (e.g., Tassoul 1980): 



Here v n \ is the cyclic frequency of a mode of order n and degree I, and a is in general a slowly 
varying function of frequency. The departures from the simple relation in Eq. (6) and Eq. (9) 
provide diagnostic information about the stellar interior. 

In the simplest approximation, the modes are uniformly spaced in radial order. The so-called 



T c (t)/R Q (t) = constant 



(5) 



3.5. Large and small p-mode separations 




(6) 
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"large separation," A, is approximately constant, 

A(n,l) ~ u ny i - v n -\,l- (7) 

The large separation A can be related to the sound speed by the formula 

R dr N 1 



In Eq. (8), the integration variable is the distance r from the center of the star; the range of 
integration is from the solar center to the surface radius R. The sound speed is denoted by c. 
According to this equation, A is the inverse of the sound travel time across a stellar diameter. 
Thus as the the sound speed decreases with increasing solar age, the sound travel time increases 
and the large separations decrease. It also follows from Eq. (6) that the frequencies of modes whose 
degrees have the same parity approximately coincide, 

v ni ~ V n -U+2 ■ (9) 
The departure from Eq. (9), the "small frequency separation," is defined by the relation 

dnl = Vnl-V n -H + 2, (10) 

and is sensitive to conditions in the stellar core. The small separation d n \ depends on the gradient 
of the solar sound speed, particularly near the stellar core. The average over n of d n i is proportional 
to a quantity D given by 

„ 1 f R dc dr . . 

D --irr- \ i — . 11 

where xq is a suitable reference value of n + l/2 (Scherrer et al. 1983; Christensen-Dalsgaard 1988). 
The average small frequency separations provide a measure of the age of the star (see, e.g., Ulrich 
1986; Christensen-Dalsgaard 1988; Guenther 1991). The solar sound speed decreases with age. 
The main decrease occurs near the core, in large part due to an increase with age of the mean 
molecular weight as a result of hydrogen being converted to helium. The sound-speed develops a 
local minimum at the center. Thus dc/dr is positive near the center of the star (in contrast to 
the negative gradient elsewhere), which reduces the value of the integral in Eq. (11). The small 
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frequency separations are more sensitive to age than the large frequency separations, although the 
large separations also decrease with increasing age. 

Table 4 lists for our standard solar model the calculated unweighted averages (denoted by 
angular brackets) of the large and small frequencies separations as a function of solar age. The 
large separations were averaged over modes with I = 0, 1, 2 and 3 and n between 10 and 22. 
This range of n was chosen because it is the observed range for most data sets. The set of small 
separations given in the third column of Table 4 (d n o) was calculated for I = and I = 2 modes; 
the small separations in the fourth column of Table 4 (d n i) were calculated for I = 1 and / = 3 
modes. The small separations were averaged over the same range of n as the large separations. 

As expected, both large and small separations decrease with time, with small separations 
decreasing much more rapidly than large separations. The decrease is slightly larger for the I = 0, 2 
separations than for the I = 1,3 separations because the I = modes are more sensitive to core 
conditions than the / = 1 modes. 

Table 4 also shows the results obtained from different helioseismology observations. References 
for the observations are given in § 6, which is where we discuss sound speeds. We do not include 
in Table 4 results for the large and small separations for LOWL1, because for this data set there 
are not enough low degree modes to form a robust average over the range of selected n. 

The calculated values for a solar age of 4.57 Gyr should be compared with the observed values 
listed in Table 4. The calculated large separation is slightly different from the observed values, 
which may reflect an uncertainty in the detailed physics of the surface layers. The uncertainty 
in the physics of the surface layers introduces errors in the calculated model frequencies. These 
frequency-dependent errors are not completely removed in the process of calculating the large 
separations from the standard solar model. The surface uncertainties are removed more successfully 
in the calculation of the small separations. 
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Table 4: The average large and small frequency separations of low degree modes of a solar model 
as a function of solar age. Also listed are the observed values for the the Sun. All splittings are 
given in units of ^Hz. 



Age (Gyr) 


(A(n,Z)> 


(dno) 


(d n i) 


u.uuu 


170.57 


19.65 


31 


.34 


U.547 


1 f* A C\ O 

154. \)6 


17.84 


28. 


.76 


1.083 


loU.zz 


16.47 


27. 


.55 


1.888 


154.87 


1 r or 

15.25 


25 


.08 


Z.Z40 


10Z.DZ 


14. /U 


24. 


.29 


2.871 


147.98 


13.42 


22 


.29 


3.408 


144.52 


12.60 


21 


.10 


4.034 


140.16 


11.67 


19 


,68 


4.570 


136.10 


10.57 


17. 


.97 


5.017 


132.87 


9.92 


16. 


.95 


5.464 


129.57 


8.98 


15. 


.53 


5.911 


126.19 


8.37 


14. 


.50 


6.358 


122.82 


7.45 


13 


.09 


6.805 


119.27 


6.90 


12. 


.07 


7.252 


115.75 


6.02 


10. 


.65 


Observations 










BiSON+LOWL 


135.33 


10.51 


17. 


.81 


GOLF 


135.12 


10.46 


17. 


.75 


GOLF2 


134.84 


10.22 


17. 


,18 


GONG 


134.86 


10.04 


18. 


.24 


MDI 


134.95 


10.14 


17. 


,64 



4. Variant and deviant solar models 

In this section, we describe eleven solar models, seven of which are slight variants on the theme 
of the standard solar model (see § 4.1) and four of which are deficient in one or more significant 
aspects of the physics used in their construction (see § 4.2). Nine of these solar models have 
been described in detail in BPB2000, where their helioseismological properties were investigated 
extensively. Hence the descriptions here will be brief. 

Since the first report (Davis, Harmer, & Hoffman 1968) that the solar neutrino event rate in the 
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chlorine experiment was less than the predicted rate (Bahcall, Bahcall, & Shaviv 1968), obtained 
by using the then standard solar model and the then standard electroweak theory, there have been 
many studies of deviant solar models that were designed to 'explain' the solar neutrino problem. 
The first non-standard solar model was proposed by Ezer & Cameron (1968), who suggested that 
the flux of 8 B neutrinos could be reduced if the central hydrogen abundance could be maintained 
near its initial value by continuous mixing in the solar core. (For a discussion of, and references 
to, the extensive early work on this problem, see Chapter 5 of Bahcall 1989.) With the advent 
of precise measurements of solar p-modes that extend deep into the solar interior, many non- 
standard models have been explicitly shown to be inconsistent with the inferred solar sound speeds 
(Bahcall, Pinsonneault, Basu, & Christensen-Dalsgaard 1997) or thep-mode frequencies (Guenther 
k Demarque 1997). 

We explore here the range of solar parameters predicted by various non-standard models, even 
those that are strongly inconsistent with helioseismological data. Our purpose is to set extreme 
limits on predicted solar parameters, such as the luminosity evolution or neutrino emission, rather 
than the traditional goal of avoiding new neutrino physics. Other authors have used non-standard 
solar models in connection with helioseismology for a variety of different applications, including, for 
example, constraining the cross section for the p — p interaction (Antia & Chitre 1999), limiting the 
amount of mass loss from the sun (Guzik & Cox 1995), and constraining the amount of anomalous 
energy transport by WIMPS (Christensen-Dalsgaard 1992). Turcotte & Christensen-Dalsgaard 
(1998a) (1998) have considered the impact of changes in heavy element abundances on the properties 
of some models; the pattern of effects they found are consistent with those obtained here. 

Table 5 and Table 6 summarize some of the important physical characteristics of the complete 
set of the nine solar models whose helioseismological properties we studied in BPB2000 plus two 
models with slightly lower values of Z/X that are studied here for the first time. We present in 
Table 5 for each of the models the central temperature, density, and pressure; the helium and heavy 
element mass fractions in the initial model; the helium and heavy element mass fractions in the 
solar center. We give in Table 6 the helium and heavy element abundances at the solar surface; the 
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mixing length parameter, a; and the radius and temperature at the base of the convective zone, 
as well as the mass included in the convective zone. All quantities are shown for the model at the 
present epoch. 

We will now define the eleven Variant and Deviant solar models. 

Table 5: Some interior characteristics of the solar models. The quantities T c (in units of 10 7 K), 
p c (10 2 gm cm -3 ), and P c (10 17 erg cm~ 3 ) are the current epoch central temperature, density, and 
pressure; Y and Z are the helium and heavy element mass fractions, where the subscript 'init' 
denotes the zero-age main sequence model, and the subscript 'c' denotes the center of the solar 
model. 



Model 




r c 


Pc 






Pc 


^mit 


Z'mit 




Y c 


Z c 


Standard 


15 


.696 


152 


.7 


2 


.342 


0.2735 


0.0188 


0. 


.6405 


0.0198 


NACRE 


15 


.665 


151 


.9 


2 


.325 


0.2739 


0.0188 


0. 


6341 


0.0197 


ASOO 


15. 


619 


152. 


.2 


2 


.340 


0.2679 


0.0187 


0. 


6341 


0.0197 


GN93 


15. 


.729 


152 


.9 


2 


.342 


0.2748 


0.02004 


0. 


.6425 


0.02110 


Pre-M.S. 


15 


.725 


152 


.7 


2 


.339 


0.2752 


0.02003 


0. 


.6420 


0.02109 


Rotation 


15 


.652 


148. 


.1 


2 


.313 


0.2723 


0.01934 


0. 


,6199 


0.02032 


Radius78 


15. 


.729 


152. 


.9 


2 


.342 


0.2748 


0.02004 


0. 


6425 


0.02110 


Radiussos 


15. 


.728 


152 


.9 


2 


,341 


0.2748 


0.02004 


0. 


6425 


0.02110 


No Diffusion 


15. 


.448 


148. 


.6 


2 


.304 


0.2656 


0.01757 


0. 


.6172 


0.01757 


Old physics 


15. 


.787 


154. 


.8 


2 


.378 


0.2779 


0.01996 


0. 


6439 


0.02102 


s u = o 


15. 


.621 


153 


.5 


2 


.417 


0.2722 


0.02012 


0. 


.6097 


0.02116 


Mixed 


15. 


.189 


90.68 


1. 


.728 


0.2898 


0.02012 


0. 


3687 


0.02047 



4.1. Variant solar models 

The NACRE model was constructed using the same input physics as our Standard model except 
that we use for the NACRE model the charged particle fusion cross sections recommended in the 
NACRE compilation (Angulo et al. 1999). We specify the fusion cross sections used for this model 
more fully in § 5.1.3. The model GN93 was considered our Standard model in BPB2000 and differs 
only in the adopted value of Z/X = 0.0245 (Grevesse & Noels 1993) from the current Standard 
model (see § 2) which has Z/X = 0.0230 (Grevesse & Sauval 1998). The model ASOO is the same as 
the two models described above except that it has lower heavy element abundance Z/X = 0.0226 
(Asplund 2000). As a consequence of a more detailed calculation of the solar atmosphere, Asplund 
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Table 6: Some characteristics of the convective zones of solar models at the current epoch. Here Y s 
and Z s are the surface helium and heavy element abundances, a is the mixing length parameter, 
-R(CZ) and T(CZ) are the radius and temperature at the base of the convective zone, and M(CZ) 
is the mass in cluded within the convective zone. 



Model 




Y a 




Z s 




a 


R(CZ) 
(Re) 


M(CZ) 
(M e ) 


T(CZ) 
(10 6 K) 


Standard 


0. 


.2437 





.01694 


2 


.04 


0.7140 


0.02415 


2.18 


NACRE 


0. 


.2443 





.01696 


2 


.04 


0.7133 


0.02451 


2.19 


ASOO 


0. 


.2386 





.01684 


2 


.05 


0.7141 


0.02394 


2.18 


GN93 


0. 


.2450 





.01805 


2 


.06 


0.7124 


0.02457 


2.20 


Pre-M.S. 


0, 


.2455 





.01805 


2 


.05 


0.7127 


0.02443 


2.20 


Rotation 


0. 


,2483 





.01797 


2 


.03 


0.7144 


0.02388 


2.15 


Radius78 


0. 


.2450 


0. 


.01806 


2 


.06 


0.7123 


0.02461 


2.20 


Radiussos 


0. 


,2450 





.01806 


2 


.06 


0.7122 


0.02467 


2.20 


No Diffusion 


0. 


.2655 





.01757 


1 


.90 


0.7261 


0.02037 


2.09 


Old physics 


0. 


,2476 





.01796 


2 


.04 


0.7115 


0.02455 


2.21 


<s 34 = 


0. 


.2422 





.01811 


2 


.03 


0.7151 


0.02309 


2.17 


Mixed 


0. 


.2535 





.01782 


1 


.85 


0.7315 


0.01757 


2.02 



(2000) suggests that all meteoritic abundances should be adjusted downward by 0.04 dex. All of 
the models described below, in § 4.2 as well as in this subsection, use the Grevesse and Noels (1993) 
composition mix with Z/X = 0.0245. 

Model Pre-M.S. is evolved from the pre-main sequence stage, but otherwise is the same as our 
Standard model. The model Rotation incorporates mixing induced by rotation and is a reasonable 
upper bound to the degree of rotational mixing which is consistent with the observed depletion of 
lithium in the sun 4 (Pinsonneault et al. 1999). The prescriptions for calculating this model are 
described in § 5 of Pinsonneault (1997) and in BPB2000. There has been considerable discussion 



4 The Rotation model discussed here differs somewhat from the rotation model analyzed in BPB2000 in that the 
metals heavier than CNO were inadvertently not mixed in the previous version of this model. The rotation profile 
computed from this model does not match precisely the best current estimates of the rotation profile in the inner 
regions of the sun. The case considered here corresponds to the maximum amount of mixing. Richard et al. (1996) 
consider solid body rotation, which corresponds to what is likely to be the minimum amount of mixing. Both of these 
models yield similar results for the effect of rotation on diffusion, which is the principal way that rotation affects solar 
neutrino fluxes. 
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recently regarding the precise value of the solar radius (cf. Antia 1998; Schou et al. 1997; Brown 
& Christensen-Dalsgaard 1998) and some discussion of the effects of the uncertainty in radius on 
the quantities inferred from the helioseismological inversions (cf. Basu 1998). We have therefore 
considered two models which were constructed with the same input physics as BP2000, but which 
have model radii which differ from the radius assumed in constructing the Standard model. Radius78 
has a radius of 695.78 Mm, which is the radius that has been determined from the frequencies of 
f-modes (cf. Antia 1998) and Radiussos has a radius of 695.508 Mm, which is the solar radius as 
determined by Brown & Christensen-Dalsgaard (1998), who used the measured duration of solar 
meridian transits during the 6 years 1981-1987 and combined these measurements with models of 
the solar limb-darkening function to estimate the value of the solar radius. 

All of these variant models are approximately as consistent with the helioseismological evidence 
as the Standard model (see BPB2000). For example, the rms sound speed differences between the 
variant models and BP2000 are: 0.03% (Pre-M.S.), 0.08% (Rotation), 0.15% (Radius 78 ), and 0.03% 
(Radiussos)- The average difference (rms) between the four variant models and the Standard model 
is 0.07% . We shall see in § 5 that the differences predicted by these models for the important 
neutrino fluxes are all less than 5%. 

4.2. Deviant solar models 

The model 'Old physics' was constructed using the old Yale equation of state (cf. Guenther et 
al. 1992), supplemented with the Debye-Hiickel correction (cf. Bahcall, Bahcall, & Shaviv 1968) 
and older OPAL radiative opacities (Iglesias, Rogers, & Wilson 1992; Kurucz 1991). The model 
includes helium and heavy element diffusion and nuclear reaction cross section data in the same 
way as our Standard model. The S34 = model was calculated assuming that the rate of the 
3 He(a,7) 7 Be reaction is zero, which implies that no 7 Be or 8 B neutrinos are produced. In the 
Standard solar model, about 12% of the terminations of the p — p chain involve the 3 He(a, 7) 7 Be 
reaction, whose rate is proportional to S34. The No Diffusion model does not include helium or 
heavy-element diffusion and therefore represents the state-of-the art in solar modeling prior to 1992 
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(cf. Bahcall & Ulrich 1988; Bahcall & Pinsonneault 1992; Proffitt 1994). The model Mixed has 
an artificially mixed core, with the inner 50% by mass (25% by radius) required to be chemically 
homogeneous at all times. This model was constructed to be similar to the prescription of Cumming 
& Haxton (1996), who changed by hand the 3 He abundance as a function of radius in the final 
BP95 (Bahcall & Pinsonneault 1995) solar model in order to minimize the discrepancy between 
measurements of the total event rates in neutrino experiments and the calculated event rates. 
Cumming and Haxton did not calculate the time evolution of their model. 

We showed in BPB2000 that the Mixed, No Diffusion, and S34 = models are strongly 
disfavored by helioseismological data. We use these deviant (or deficient) models here to test 
the robustness of the discrepancies between solar model predictions and solar neutrino measure- 
ments. (For references and discussions to earlier work on these deviant models, see Bahcall 1989, 
Christensen-Dalsgaard 1995, Bahcall, Pinsonneault, Basu, Sz Christensen-Dalsgaard 1997,, Guen- 
ther & Demarque 1997, and Cox, Kidman, and Newman 1985). We have already seen in Figure 2 
that the luminosity evolution predicted by solar models is essentially the same for the deviant 
models and for the Standard solar model. 



5. Neutrino Physics 

This section presents neutrino fluxes and predicted event rates in different solar neutrino 
experiments; the fluxes and event rates are calculated using the eleven solar models described in 
§ 4. We also give other aspects of the solar model that influence the interpretation of solar neutrino 
experiments. 

As described in § 2.2, we present at http://www.sns.ias.edu/~jnb a detailed numerical table 
which gives the fraction of each of the eight important neutrino fluxes that is produced in each 
spherical shell. These neutrino production fractions are important for calculating the effect of MSW 
(matter) oscillations in the sun, but will not be discussed further here. 

Section 5.1 presents the fluxes of electron type neutrinos that are produced in the sun according 
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to the Standard solar model and to the eight variant solar models considered in this paper. We 
compare the results of these predictions with measurements of the total rates in the chlorine solar 
neutrino experiment, the SAGE and GALLEX + GNO gallium experiments, and the Kamiokandc 
and Super-Kamiokande water Cherenkov detectors. We also describe in this section how we estimate 
the errors (asymmetric in some cases) on the predicted fluxes and event rates. 

In § 5.2, we present detailed results for the electron number density versus solar radius. The 
MSW conversion of electron type neutrinos to other neutrino types depends upon the radial distri- 
bution of the electron number density. In previous presentations of the Standard model, we have 
not given the electron number density in the outer parts of the sun with sufficient precision to 
permit calculations of a subset of the currently allowed matter transitions. Section 5.4 presents the 
calculated time-dependences of the most important solar neutrino fluxes, the pp, 7 Be, 8 B, and 13 N 
neutrino fluxes. 

5.1. Neutrino fluxes and experimental event rates at the current epoch 

We present in § 5.1.1 the neutrino fluxes and experimental event rates predicted by the Stan- 
dard model and contrast these results with the observed rates. In § 5.1.2, we describe the procedures 
and the ingredients used to calculate the uncertainties in the neutrino fluxes and the event rates. 
We compare in § 5.1.4 the calculated neutrino fluxes and experimental event rates for the eight 
variant and deviant solar models with the results for the Standard model and with the measured 
solar neutrino event rates. 

5.1.1. Standard model 

Table 7 gives the neutrino fluxes and their uncertainties for our Standard solar model. In 
order to obtain the most precise values that we can for the predicted fluxes, we have recomputed 
the Standard model discussed elsewhere in this paper. We use in this subsection the most recently 
determined absolute value for the solar luminosity, 3.842 x 10 33 ergs -1 (Frohlich & Lean 1998, 
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Table 7: Standard Model Predictions (BP2000): solar neutrino fluxes and neutrino capture rates, 
with la uncertainties from all sources (combined quadratically) . The tablulated fluxes correspond 
to a present-day solar luminosity of 3.842 x 10 33 ergs -1 . The observed capture rates are: 2.56 ±0.23 
SNU [chlorine (Lande 2000)] and 74.7±5.0 SNU [combined SAGE and GALLEX plus GNO (Hampel 
et al. 1999; Abdurashitov et al. 1999; Bellotti 2000)]. The 8 B flux measured by the Super- 
Kamiokande experiment is 2.40 ± 0.03(stat) ^^(syst.) cnrV 1 (Suzuki 2000). 



Source 


Flux 


CI 


Ga 


Li 




(10 10 cm^s" 1 ) 


(SNU) 


(SNU) 


(SNU) 


PP 


5.95(1.001851) 


0.0 


69.7 


0.0 


pep 


i.4o x io- 2 (i.ooi8:8^) 


0.22 


2.8 


9.2 


hep 


9.3 x 10~ 7 


0.04 


0.1 


0.1 


7 Be 


4.77 x IO -1 (1.0018 J8) 


1.15 


34.2 


9.1 


8 B 


5.05 x io- 4 (i.ool8:?8) 


5.76 


12.1 


19.7 


13 N 


5.48 x IO -2 (I.00l8:i7) 


0.09 


3.4 


2.3 


15Q 


4.80 x io~ 2 (i.ool8:?g) 


0.33 


5.5 


11.8 


17 F 


5.63 x 10~ 4 (1.0018:25) 


0.0 


0.1 


0.1 


Total 




7.6I 1 : 3 


1281? 


co q+6.5 
oz -°-6.0 



Note. — The cross sections for neutrino absorption on chlorine are from Bahcall et al. (1996); the cross sections 
for gallium are from Bahcall (1997); the cross sections for 7 Li are from Bahcall (1989) and Bahcall (1994). 

Crommelynck, Fichot, Domingo, & Lee 1996), which is 0.2% smaller than the value used in the 
model calculations discussed elsewhere in this paper. The largest changes due to the 0.2% decrease 
in the solar luminosity are a 2% decrease in the 8 B neutrino flux and a 1% decrease in the 7 Be 
neutrino flux. All other quantities calculated in this paper are changed by negligible amounts. 

The adopted uncertainties in different input parameters are given in Table 2 of Bahcall, Basu, 
& Pinsonneault (1998), which we refer to hereafter as BP98. We also present in Table 7 the 
calculated event rates in the chlorine, gallium, and lithium experiments. The rate in the Super- 
Kamiokande experiment is usually quoted as a fraction of the best-estimate theoretical flux of 8 B 
neutrinos, assuming an undistorted (standard) energy spectrum. 

Table 8 compares the predictions of the combined standard model, i. e., the standard solar 
model (BP2000) and the standard electroweak theory (no neutrino oscillations), with the results 
of the chlorine, GALLEX + GNO, SAGE, Kamiokande, and Super-Kamiokande solar neutrino 
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Table 8: Solar Neutrino Rates: Theory versus Experiment. The units are SNU (1CT 36 interactions 
per atom per sec) for the radiochemical experiments: Chlorine, GALLEX + GNO, and SAGE. 
The unit for the 8 B and hep fluxes are, respectively, 10 6 cm -2 s _1 and 10 3 cm~ 2 s _1 . The errors 
quoted for Measured/BP2000 are the quadratically combined uncertainties for both BP2000 and 
the Measured rates. For simplicity in presentation, asymmetric errors were averaged. References to 
the experimental results are given in the text and in Lande (2000), Bellotti (2000), Gavrin (2000), 
Fukuda et al. (1996), and Suzuki (2000) for the chlorine, GALLEX + GNO, SAGE, Kamiokande, 
and Super-Kamiokande results. 



Experiment 


BP2000 


Measured 


Measured/BP2000 


Chlorine 


7.6 + {i 


2.56 ±0.23 


0.34 + 0.06 


GALLEX + GNO 


128+ 7 


74 1+ 6 - 7 


0.58 + 0.07 


SAGE 


128+ 7 


75.4+^; 8 


0.59 + 0.07 


8 B-Kamiokande 


5.05 [l.OO+l™] 


2.80 [1.00 ±0.14] 


0.55 + 0.13 


8 B-Super-Kamiokande 


5.05 [l.00+18:?g] 


2.40 [l.OO+18;8t] 


0.48 ± 0.09 


hep-Super-Kamiokande 


9.3 


11.3(1+0.8) 


~ 1 



experiments. The observed rate in the chlorine experiment is 2.56 ± 0.23 SNU (Lande 2000; Davis 
1994; Cleveland et al. 1998), which is to be compared to the calculated value of 7.61J' 3 SNU. This 
discrepancy between calculated and observed neutrino capture rates has been approximately the 
same for more than three decades (cf. Bahcall, Bahcall, & Shaviv 1968; Davis, Harmer, & Hoffman 
1968; Bahcall 1989). 

The average of the SAGE (Abdurashitov et al. 1999; Gavrin 2000) and the GALLEX (Hampel 
et al. 1999) plus GNO (Bellotti 2000) results is 74.7 ± 5.0 SNU, which is more than 6cr away from 
the calculated standard rate of 1281? SNU. 

After 1117 days of data acquisition, the flux of 8 B neutrinos measured by Super-Kamiokande is 
[2.40 ± 0.03(stat) 18;87(syst.)]10~ 6 chi^V 1 (Suzuki 2000; Fukuda et al. 1998), which corresponds 
to 0.475 of the BP2000 predicted flux. 

Comparing the second and third columns of Table 8, we see that the predictions of the com- 
bined standard model differ by many standard deviations from the results of the solar neutrino 
experiments. 

The flux of hep neutrinos was calculated for BP2000 using the most recent theoretical eval- 
uation by Marcucci et al. (2000a,b) of the cross section factor So (hep), which is 4.4 times larger 
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than the previous best estimate. The most recent preliminary report (after 1117 days of data) of 
the Super-Kamiokande collaboration (Suzuki 2000) is that the hep flux observed in their v e — e 
scattering experiment is 5.4 ± 4.5 times the best-estimate from BP98. Since the BP2000 estimate 
of the hep flux is a factor of 4.4 times larger than the flux quoted in BP98, the best-estimate the- 
oretical hep flux now agrees with the best-estimate experimental hep flux measurement, although 
we do not attach much significance to this agreement since we cannot quote an uncertainty on the 
theoretical estimate (see discussion of the hep reaction in § 5.1.2). 

The event rates predicted by BP2000 for the chlorine, gallium, and Super-Kamiokande solar 
neutrino experiments are within two percent of the rates predicted for the BP98 standard solar 
models. As far as these experiments are concerned, the effects of using the improved heavy element 
composition essentially cancels the effect of correcting the error in the opacity interpolation (see 
§2.1). The difference in the 7 Be flux predicted by BP98 and BP2000 is only 0.6%. The 7 Be flux 
will be measured by the BOREXINO solar neutrino experiment. 



5.1.2. Calculated uncertainties 

We have calculated the uncertainties in the neutrino fluxes and in the experimental event 
rates by including the published errors in all experimental quantities and by taking account of 
the correlated influence of different input parameters using the results of detailed solar model 
calculations. The procedure for calculating the uncertainties has been developed over the past three 
decades and is described in detail in Chapter 7 of Bahcall (1989) (see also Bahcall & Pinsonneault 
1992, 1995, and Bahcall, Basu, & Pinsonneault 1998). 

In order that the reader can see the specific implementation of the uncertainty calculations, 
we are making available the exportable fortran code that evaluates the rates in different neutrino 
experiments and also calculates the uncertainties in the individual neutrino fluxes and experimental 
rates. The code, exportrates.f, is available at http://www.sns.ias.edu/~jnb/SNdata/sndata.html . 

The uncertainties in the nuclear fusion cross sections (except for hep, see below) were taken 
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from Adelberger et al. (1998), the neutrino cross sections and their uncertainties are from Bahcall 
(1994, 1997) and Bahcall et al. (1996), the luminosity and age uncertainties were adopted from 
Bahcall & Pinsonneault (1995), the la fractional uncertainty in the diffusion rate was taken to be 
15% (Thoul, Bahcall, & Loeb 1994) and the opacity uncertainty was determined by comparing the 
results of fluxes computed using the older Los Alamos opacities with fluxes computed using the 
modern Livermore opacities (Bahcall & Pinsonneault 1992). 

We follow the discussion in Bahcall & Pinsonneault (1995) and adopt a la uncertainty in the 
heavy element abundance of 

a{Z/X) = ±0.061 x (Z/X). (12) 

This uncertainty spans the range of values recommended by Grevesse (1984), Grevesse & Noels 
(1993), and Grevesse Sz Sauval (1998) over the fourteen year period covered by the cited Grevesse 
et al. review articles. The uncertainty adopted here is about twice as large as the uncertainty 
recommended by Basu & Antia (1997) based on their helioseismological analysis. In support of the 
larger uncertainty used in this paper, we note that the difference between the Grevesse & Noels 
(1993) values of Z/X = 0.0230 is la according to Eq. (12). 

We include for the first time in this series of papers the uncertainty in the small 17 F neutrino 
flux due to the uncertainty in the measured So-factor for the reaction 16 0(p, 7) 17 F. We use the la 
uncertainty 18.1% estimated by Adelberger et al. (1998). It was an oversight not to include this 
uncertainty in our previous calculations. 

The only flux for which we do not quote an estimated uncertainty is the hep flux (see Table 7) . 
The difficulty of calculating from first principles the nuclear cross section factor So (hep) is what 
has caused us not to quote in this series of papers an uncertainty in the hep flux (see discussion in 
Bahcall 1989 and in Bahcall & Krastev 1998). The hep reaction is uniquely difficult to calculate 
among the light element fusion reactions since the one-body and two-body contributions to the 
reaction rate are comparable in magnitude but opposite in sign, so that the net result is sensitive 
to a delicate cancellation. Also, two-body axial currents from excitations of A isobars are model 
dependent. In addition, the calculated rate is sensitive to small components in the wave function, 
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particularly D-state admixtures generated by tensor interactions. These complications have been 
discussed most recently and most thoroughly by Marcucci et al. (2000b). 

The calculated errors are asymmetric in many cases. These asymmetries in the uncertainties 
in the neutrino fluxes and experimental event rates result from asymmetries in the uncertainties 
of some of the input parameters, for example, the important pp, 7 Be + p, and 14 N + p fusion 
reactions and the effect of excited states on neutrino absorption cross sections. To include the 
effects of asymmetric errors, the code exportrates.f was run with different input representative 
uncertainties and the different higher (lower) rates were averaged to obtain the quoted upper 
(lower) limit uncertainties. 

5.1.3. NA CRE charged particle fusion rates 

In order to estimate the systematic uncertainties associated with different treatments of the 
nuclear fusion reactions, we have constructed a solar model that is the same as the Standard Model 
discussed in § 5.1.1, except that we have used the charged particle fusion cross sections recommended 
in the NACRE compilation (Angulo et al. 1999) rather than the fusion cross sections determined 
by Adelberger et al. (1998). We will refer to this solar model as the NACRE model. 

The low energy cross section factors, So, that are recommended by the NACRE collaboration 
and by Adelberger et al. agree within their stated la uncertainties for all of the fusion reactions that 
are important for constructing a solar model. The only important solar nuclear reactions for which 
the NACRE collaboration did not recommend interaction rates are the electron capture reactions 
that produce the 7 Be and the pep neutrinos; the NACRE collaboration also did not provide energy 
derivatives for the cross section factors of the CNO reactions. Wherever the data necessary for 
computing solar fusion rates was not available in the NACRE compilation, we continued to use the 
Adelberger et al. (1998) recommended values in computing the NACRE model. 

Table 9 gives the calculated neutrino fluxes and capture rates predicted by the NACRE solar 
model. In all cases, the fluxes for the NACRE solar model agree with the fluxes calculated with 
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Table 9: Neutrinos with NACRE reaction rates. The solar neutrino fluxes and neutrino capture 
rates that were calculated with the NACRE fusion cross sections (Angulo et al. 1999) are given in 
the table. The only difference between the model used in this calculation and the Standard Model, 
whose fluxes are given in Table 7, is that the Adelberger et al. (1998) fusion cross sections were 
replaced by the NACRE cross sections for all reactions for which the NACRE collaboration quoted 
zero-energy cross section factors, Sq. 



Source 


Flux 


CI 


Ga 


Li 




(10 10 cm-V 1 ) 


(SNU) 


(SNU) 


(SNU) 


PP 


5.96(1.00^}) 


0.0 


69.8 


0.0 


pep 


1.39 x 10- 2 (l.OOlffi) 


0.22 


2.8 


9.1 


hep 


9.4 x 10~ 7 


0.04 


0.1 


0.1 


7 Be 


4.81 x lO" 1 (l.OOlo 


1.15 


34.5 


9.1 


8 B 


5.44 x lO" 4 (1.001°;?°) 
4.87 x 10~ 2 (l.OOlg;^) 
4.18 x lO" 2 (l.OOlJujj) 
5.30 x 10~ 4 (l.00±g;||) 


6.20 


13.1 


21.2 


13 N 


0.08 


2.9 


2.1 


15Q 


0.28 


4.7 


10.3 


17p 


0.0 


0.1 


0.1 


Total 




8.01}; 4 


1281? 


52-0±g;| 



Note. — The cross sections for neutrino absorption on chlorine are from Bahcall et al. (1996); the cross sections 
for gallium are from Bahcall (1997); the cross sections for 7 Li are from Bahcall (1989) and Bahcall (1994). 

the Standard solar model to well within the la uncertainties in the Standard Model fluxes. The 
7 Be flux from the NACRE model is 1% larger than for the Standard model and the 8 B flux is 
8% higher. The chlorine capture rate predicted by the NACRE model is 5% higher than for the 
Standard model; the predicted rates for the NACRE model and the Standard model differ by less 
than 1% for the gallium and lithium experiments. 

We conclude that this estimate of the systematic uncertainties due to the relative weights given 
to different determinations of nuclear fusion cross sections suggests likely errors from this source 
that are significantly smaller than our quoted la errors for the Standard model neutrino fluxes 
and event rates (cf. Table 7). Similar conclusions have been reached by Morel, Pichon, Provost, & 
Berthomieu (1999) in an independent investigation (see also Castellani et al. 1997). 
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Table 10: Neutrino fluxes from twelve solar models. The table lists the calculated neutrino fluxes 
and event rates for the seven variant solar models that are discussed in § 4.1 and for the four 
deficient solar models that are discussed in § 4.2, and compares the results with the Standard 
model fluxes. The models in rows two through eight are all variants on the Standard model (first 
row). The last four models are deficient in some important aspect of the physics used in their 
construction and do not provide good fits to the helioseismological data. 



Model 


PP 


pep 


hep 


Be 


8r 

D 


13 AT 
IN 


i5n 
u 


17T? 
r 


L/l 


(jSl 








[Hid) 




[hit)) 










/CWTT\ 


Standard 


5.96 


1.40 


9.3 


4.82 


5.15 


5.56 


4.88 


5.73 


7.7 


129 


NACRE 


5.97 


1.39 


9.4 


4.85 


5.54 


4.93 


4.24 


5.39 


8.1 


129 


AS00 


5.99 


1.41 


9.4 


4.62 


4.70 


5.25 


4.56 


5.33 


7.1 


126 


GN 93 


5.94 


1.39 


9.2 


4.88 


5.31 


6.18 


5.45 


6.50 


8.0 


130 


Pre-M.S. 


5.95 


1.39 


9.2 


4.87 


5.29 


6.16 


5.43 


6.47 


7.9 


130 


Rotation 


5.98 


1.40 


9.2 


4.68 


4.91 


5.57 


4.87 


5.79 


7.4 


127 


Radius78 


5.94 


1.39 


9.2 


4.88 


5.31 


6.18 


5.45 


6.50 


8.0 


130 


Radiussos 


5.94 


1.39 


9.2 


4.88 


5.31 


6.18 


5.45 


6.50 


8.0 


130 


No Diffusion 


6.05 


1.43 


9.6 


4.21 


3.87 


4.09 


3.46 


4.05 


6.0 


120 


Old physics 


5.95 


1.41 


9.2 


4.91 


5.15 


5.77 


5.03 


5.92 


7.8 


130 


S 3 a = 


6.40 


1.55 


10.1 


0.00 


0.00 


6.47 


5.64 


6.70 


0.8 


89 


Mixed 


6.13 


1.27 


6.2 


3.57 


4.13 


3.04 


3.05 


3.61 


6.1 


115 



5.1.4- Variant and deviant models 



Table 10 compares the calculated neutrino fluxes for the seven variant solar models described 
in § 4.1 and for the four deficient solar models described in § 4.2 with the fluxes obtained for the 
Standard solar model. The range of the 8 B neutrino fluxes among the seven standarddike models 
(rows 1-7 of Table 10) is only ±7%, a factor of two or three smaller than the uncertainty (due to 
other sources) in the calculated 8 B flux of the Standard model (see Table 7). The spread in 7 Be 
flux is ±3%; the range in 37 C1 and 71 Ga rates is ±0.45 SNU and ±2 SNU, respectively. 

The deviant models listed in Table 10 are all deficient in some important aspect of the physics 
used in their calculation. The No Diffusion, Old Physics, S34 = 0, and Mixed models all give such 
bad agreement with helioseismology (compared to the first seven models) that the comparison with 
the data cannot be made on the same scale as for the standarddike models (see BPB2000). The 
Old Physics model gives a rms difference between the calculated and measured sound speeds that 
is more than twice as large as when the Standard model is used. For the 534 = model, the 
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rms discrepancy is seven times worse in the solar interior than for the standard model; for the No 
Diffusion model, the disagreement is about seven times worse averaged over the sun than for the 
standard model. The Mixed model is the worst of all; the rms disagreement in the solar core is 
about 25 times larger than for the standard solar model. 

Even if one is willing to consider solar models that predict sound speeds which deviate so 
drastically from the measured helioseismological values, the four deviant solar models do not de- 
scribe well the solar neutrino data. For the Old Physics, No Diffusion, and Mixed solar models, 
the predicted rates for the gallium experiments lie between 115 SNU and 130 SNU, which is to 
be contrasted with the observed rate of 75 ± 5 SNU, which is at least an 8a discrepancy. For the 
S34 = model, the predicted rate in the chlorine experiment is 0.79 SNU, which is also about eight 
standard deviations less than the observed value, 2.56 ± 0.23 SNU. 



The probability of converting an electron type neutrino to a muon or tau neutrino in the 
Sun depends upon the profile of the electron number density as a function of solar radius. For 
particular values of the electron density, neutrino energy, and neutrino mass, neutrinos can be 
resonantly converted from one type of neutrino to another. The Mikheyev-Smirnov resonance 
occurs if the electron density at a radius r satisfies 



where n e is the electron number density measured in cm~ 3 , Na is Avogadro's number, 6y is the 
neutrino mixing angle in vacuum, |Am 2 | is the absolute value of the difference in neutrino masses 
between two species that are mixing by neutrino oscillations, and E is the neutrino energy. 

Figure 8 and Table 11 give the electron number density as a function of solar radius for the 
Standard solar model (BP2000). A much more extensive numerical file of the electron number 
density versus radius is available at http://www.sns.ias.edu/~jnb; this file contains the computed 
values of the electron number density at 2493 radial shells. 



5.2. The electron number density 




(13) 
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log(n e /N A ) vs. R/R 
BP2000 




Fig. 8. — The electron number density, n e , versus solar radius for the Standard solar model 
(BP2000). The straight-line fit shown in Figure 8 is an approximation, Eq. (14), given by Bahcall 
(1989). Equation (14) has been used previously in many analyses of matter effects on solar neutrino 
propagation. Precise numerical values for n e are available at 



Table 11: The electron number density versus radius for the standard solar model. The tabulated 
values are log(n e /iV^), where n e is measured in number per cm 3 and Na is Avogadro's number. 
A more extensive numerical file of n e , contain electron number densities at 2493 radial shells, is 
available at http://www.sns.ias.edu/~jnb. 
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We see from Figure 8 that for typical values of the neutrino parameters that allow the so-called 
LMA and SMA MSW solutions which fit all of the currently available solar neutrino data (e.g., 
Bahcall, Krastev, h Smirnov 1998), the experimentally most important 8 B neutrinos (E > 5 MeV) 
satisfy the resonance condition, Eq. (13), at radii that are smaller than the radius of the convective 
zone. For the so-called LOW MSW solutions and for all MSW solutions with Oy ~ 7r/4, the 
resonance radius falls in the outer part of the sun. 

We have not previously published accurate values for the electron density in the outer parts 
of the sun, r > O.8i? . The straight line in Figure 8 is an approximation to the electron 
number density in the standard solar model of Bahcall & Ulrich (1988) (see Bahcall 1989). This 
approximation, which has been used for over a decade by different groups analyzing solar neutrino 
data, is 

n e /N A = 245exp(-lO.54i?/ii ) cm" 3 . (14) 

Figure 8 shows that the approximation given in Eq. (14) fails badly in the outer regions of the 
sun. Recently, several different analyses have been published using the electron number den- 
sity shown in Figure 8 (or, more precisely, the computer file for n e (r) which is available at 
http://www.sns.ias.edu/~jnb ). 

5.3. The number density of scatterers of sterile neutrinos 

The effective density of particles for interacting with sterile (right handed) neutrinos is not the 
electron number density discussed in the previous subsection, but rather n ste riie (see Mikheyev Sz 
Smirnov 1986; Lim & Marciano 1988; Barger et al. 1991), where 

^•sterile = 0.5 X Tlneutronsj (1^) 

and n neu trons is the number density of neutrons. Since nearly all of the neutrons in the sun are 
either in 4 He or in heavier elements with Z ~ A/2, it is easy to derive an analytic expression that 



relates n ster ii e and n e . One obtains 



nstenle = U e (^f)) ■ (16) 

Figure 9 and Table 12 give the radial distribution of ra s t er iie in. the Standard solar model. The 
functional form of re ster ii e (?") is similar to the function form of n e {r). The straight line in Figure 9 is 
given by an equation of the same form as Eq. (14) that describes n e (r) except that the coefficient 

for Tlgterile 

/N A is 223 (instead of 245 for n e /N A ). 

The number density n s t er iie is about 25% smaller than the electron number density, n e , in the 
center of the sun, where helium is most abundant. In the central and outer regions of the sun, 
^sterile is about 9% less than n e . Since the slopes of the straight-line fitting functions are the same 
in Figure 8 and Figure 9, the neutrino survival probabilities are the same for sterile and for active 




Fig. 9. — The number density, n ster iie, of scatterers of sterile neutrinos versus solar radius for the 
Standard solar model (BP2000). The straight line in Figure 9 is given by an equation of the same 
form as Eq. (14) except that the coefficient for 7i s teriie 

/N A is 223 (instead of 245 for n e /N A ). 
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Table 12: The sterile number density, n steri i e , versus radius for the standard solar model. The 
tabulated values are log(n ste riie,MOl) 5 where n ster iie is measured in number per cm 3 and Na is 
Avogadro's number. A more extensive numerical file of n ster iie containing values of ra s teriie at 2499 
radial shells, is available at http://www.sns.ias.edu/~jnb. 
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neutrinos as long as the adiabatic approximation is valid (see, e.g., § 9.2 of Bahcall 1989). 



5.4. Neutrino fluxes as a function of solar age 

Figure 10 shows the most important solar neutrino fluxes, the pp, 7 Be, 8 B, and 13 N fluxes, as 
a function of solar age. The fluxes displayed in the figure were computed using the Standard model 
and are normalized by dividing each flux by its value at the present epoch, labeled by the word 
'today' in the figure (cf. Fig. 10 of Guenther and Demarque 1991 for a similar figure plotted on a 
logarithmic scale). 

The pp flux is relatively constant over the entire 8 billion years show in Figure 10. At the 
beginning of its lifetime, the pp flux is about 75% of its current value and reaches 90% of its present 
value after 2.6 billion years. At the current epoch, the flux is changing very slowly, about 4% per 
billion years. The pp flux reaches a maximum, 4% larger than its current value, at a solar age of 
6.0 billion years and then declines slowly and steadily to 96% of its present value at an age of 8 
billion years. 
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Fig. 10. — The pp, 7 Be, 8 B, and 13 N neutrino fluxes as a function of solar age. The figure shows 
the Standard model ratios of the fluxes divided by their values at 4.57 x 10 9 yr. The pp flux is 
represented by a solid line, the 7 Be flux by a line of long dashes, the 8 B flux by short dashes, and 
the 13 N flux by a dotted line. 

The 7 Be and 8 B neutrino fluxes increase monotonically and by larger amounts than the pp 
flux. Both the 7 Be and the 8 B fluxes begin with very low fluxes relative to their current values, 
14% and 3%, respectively, of their intensities at 4.57 x 10 7 yr. At a solar age of 8 billion years, the 
7 Be neutrino flux is 2.6 times larger than it is today and the 8 B neutrino flux is 7.1 times larger 
than today. At the current epoch, the 7 Be flux is increasing by about 65% per billion years and 
the 8 B flux is increasing faster, about 120% per billion years. 

The 13 N neutrino flux has the most interesting time dependence. In the first 10 8 y on the 
main sequence, the 13 N flux is much larger than its current value because 12 C has not yet been 
burned to the equilibrium value appropriate for the CNO cycle. The reaction 12 C(p, 7) 13 N occurs 
relatively often in this early stage of solar evolution and the neutrino flux from 13 N beta-decay has 
a peak value of about 11 times its current flux. The minimum 13 N flux, 33% of its present value, is 
attained at a solar age of 1.8 billion years. Thereafter, the 13 N flux increases steadily as the central 
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temperature of the solar model increases and reaches an intensity of 18 times its current value at 
a solar age of 8 billion years. 

6. Sound speeds 

Section 6.1 presents a panoramic view of the predicted Standard model sound speeds and 
compares the observations and the calculations on a scale that is relevant for interpreting solar 
neutrino experiments. In § 6.2, we compare on a zoomed- in scale the Standard model calculations 
with the results of six helioseismological measurements. The zoomed-in scale used in this subsection 
can reveal fractional discrepancies between calculations and observations that are smaller than 0.1%. 
In § 6.3, we consider a particular solar model that includes rotation in a plausible way. This model 
smoothes the composition discontinuity at the base of the convective zone, which locally improves 
the agreement with the measured sound speeds but worsens the overall rms agreement (within the 
uncertainties allowed by input data). For neutrino emission, the predictions of the rotation model 
are not significantly different from the Standard model predictions. 

Some authors (see, e.g., Guenther & Demarque 1997) compare their solar models directly with 
the p-mode oscillation frequencies rather than with inverted quantities such as the sound speed. The 
reader is referred to the Guenther and Demarque paper for a discussion of the direct comparison 
method, its application, and additional references. 

We have chosen to use the sound speed profile because the inversion process that produces 
the inferred sound speeds allows one to remove the uncertainties, common to all p-mode oscillation 
frequencies, that arise from the near-surface regions of the sun. These common uncertainties are 
due to the treatment of convection, turbulence, and non-adiabatic effects. Inversion techniques 
are designed to minimize the effects of these outer layers (see, e.g., Basu, Christensen-Dalsgaard, 
Hernandez, and Thompson 1996) Moreover, the sound speed profile summarizes in a robust way 
the results obtained for many thousands of oscillation frequencies. Finally, the inversion procedure 
allows one to isolate different regions of the sun, which is important in the context of discussions 
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regarding solar neutrinos. The neutrinos are produced deep in the solar interior. In BPB2000, we 
have discussed in detail the systematic uncertainties and assumptions related to the inversion for 
the sound speeds and for the less-accurately determined density. 



6.1. Sounds speeds: panoramic view 

Figure 11 shows the fractional differences between the calculated sound speeds for the Standard 
model and what may be the most accurate available sound speeds measured by helioseismology, the 
LOWL1 + BiSON Measurements presented in Basu et al. (1997). These sound speeds are derived 
from a combination of the data obtained by the Birmingham Solar Oscillation Network (BiSON; 

— , — , — i — i — i — , — | — , — , — , — 

Bahcall-Pinsonneault 2000 J 
LOWL1 + BiSON Measurement - 
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Fig. 11. — Predicted versus measured sound speeds. The figure shows the excellent agreement 
between the calculated sound speeds for the Standard solar model (BP2000) and the helioseismo- 
logically measured (Sun) sound speeds. The horizontal line at 0.0 represents the hypothetical case 
in which the calculated sound speeds and the measured sound speeds agree exactly everywhere in 
the sun. The rms fractional difference between the calculated and the measured sound speeds is 
0.10% for all solar radii between between 0.05i? Q and 0.95R & and is 0.08% for the deep interior 
region, r < O.25i?0, in which neutrinos are produced. 



U. 1 

0.08 h 
0.06 
0.04 - 
0.02 - 




c 

3 

m 

\ 

c 

3 

m 

l 

« -0.02 
o 

S -0.04 - 
-0.06 - 
-0.08 
-0.1 



^ 1 



7 Be lowered (lcroff ( 



-47- 



cf. Chaplin et al. 1996) and the Low-£ instrument (LOWL; cf. Tomczyk et al. 1995a, b). 



The rms fractional difference between the calculated and the measured sound speeds is 10.4 x 
10~ 4 over the entire region in which the sound speeds are well measured, O.O5i?0 < r < O.95i?0. 
In the solar core, 0.05-R© < r < 0.25i? Q (in which about 95% of the solar energy and neutrino 
flux are produced in a standard solar model), the rms fractional difference between measured and 
calculated sound speeds is 6.3 x 10 -4 . The Standard model sound speeds agree with the measured 



sound speeds to 0.1% whether or not one limits the comparison to the solar interior or averages 
over the entire sun. Systematic uncertainties ~ 3 x 10~ 4 are contributed to the sound speed profile 



by each of three sources: the assumed reference model, the width of the inversion kernel, and the 



measurement errors (see BPB2000). 



The vertical scale of Figure 11 was chosen so as to include the arrow marked " Be lowered (la off Ga). 
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Fig. 12. — The solar sound speed versus the solar radius. The figure shows the calculated solar sound 
speed versus radius for the standard solar model, Bahcall-Pinsonneault (2000). To an accuracy of 
about 0.5 km/s the calculated and the observed sound speeds are the same (see Figure 13 and 
Figure 14). 
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This arrow indicates the typical difference between solar model speeds and helioseismological mea- 
surements that would be expected if the discrepancy between the gallium solar neutrino measure- 
ments and the predictions in Table 7 were due to errors in the solar physics of the standard solar 
model (see discussion in BP98). 

Figure 12 and Table 13 give the sound speeds versus the solar radius that are calculated using 
the Standard solar model. The sound speed declines from about 500 km/s in the solar core to 
about 100 km/s at O.95i? . An extensive table of the Standard model sound speeds is available at 
http://www.sns.ias.edu/~jnb . 

Table 13: The sound speed as a function of radius for the Standard solar model. The Sound speeds 
are given in units of 100 km/s. 
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6.2. Sounds speeds: zoom in 

Figure 13 compares the results of six precise observational determinations of the solar sound 
speed with the results of our Standard solar model. The vertical scale has been expanded by a 
factor of 21 relative to Figure 11 in order to show the small but robust discrepancies between the 
calculations and the observations and to indicate the size of the differences between the various 
measurements. In the deep solar solar interior where neutrinos are produced, R/R Q < 0.25, the 
differences between the various observational determinations of the sound speed are comparable to 



-49- 



0.004 



0.002 - 



c 

3 

m 



c 

3 

m 

T3 
O 



-0.002 - 



-0.004 - 



-i — i — i — | — , — , — , — | — , — , — , — r 
Bahcall-Pinsonneault 2000 



-i — i — | — i — i r- 

L0WL1 + BiSON 

GOLF 

G0LF2 

GONG 

LOWLl 

MDI 




0.4 0.6 
R/R 



Fig. 13. — Six precise helioseismological measurements versus BP2000. The figure compares the 
fractional difference between the sound speeds calculated for the Standard solar model (BP2000) 
and the sound speeds in six helioseismological experiments. The references to the helioseismological 
data are given in the text. Systematic uncertainties due to the assumed reference model and the 
width of the inversion kernel are each ~ 0.0003 (see BPB2000). 



the differences between BP2000 and any one of the measured sets of sound speeds. 

The p— mode frequencies used in deriving the observed sound speeds shown in Figure 13 were 
obtained from a number of different sources. In addition to the LOWLl + BiSON data described 
in § 6.1, we have used data from a number of other sources. (1) Data from the first year of LOWL 
observations. The sound speed inversions are described in Basu et al. (1997) and are referred to as 
LOWLl in this paper. (2) Data from the Michelson Doppler Imager (MDI) instrument on board the 
Solar and Heliospheric Observatory (SOHO) during the first 144 days of its operation (cf. Rhodes 
et al. 1997). The results of the sound speed inversions using these data are from Basu (1998). 

(3) The frequencies obtained from the data obtained by the Global Oscillation Network Group 
(GONG) between months 4-14 of its observations. The solar sound-speeds are from Basu (1998). 

(4) Initial observation taken by the Global Oscillations at Low Frequencies (GOLF) instrument 
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on board SOHO, combined with intermediate data from MDI. We have labeled the sound speeds 
obtained as GOLF1 and the solar sound speed results can be found in Turck Chieze et al. (1997). 
(5) More recent data from GOLF (Thiery et al. 2000), combined with intermediate degree data 
obtained from the first 360 days observations by the MDI instrument (Schou et al. 1997). The 
sound-speed results are described in Basu et al. (2000). These results have been labeled as GOLF2. 

There are other helioseismological data sets that could have been used. Of these, the data by 
Toutain et al. (1998) may be the most relevant. This data set has been discussed extensively in 
Basu et al. (2000). Toutain et al. (1998) were the first to show the effect of line-profile asymmetry 
on low-degree modes. The sound-speed inversions in their paper were, however, very different 
from what had been obtained previously. Basu et al. (2000) showed that this difference results 
primarily from the two modes, I = 2, n = 6 and I = 2 n = 7, whose frequencies (and errors 
on the frequencies) were suspect. In a later paper, Bertello et al. (2000) confirmed most of the 
low frequency modes in the Toutain et al. (1998) data set, but not the two questionable modes. 
Additionally, they determined the frequencies of some very low-frequency low-degree modes. The 
sound-speed inversion of Bertello et al. (2000) is similar to the inversion obtained in Basu et al. 
(2000) and in the present paper; the errors in the inversions are also similar. The similarity in the 
errors may appear somewhat surprising at first glance given that the Bertello et al. frequencies 
have very low errors, but the similarity in the overall uncertainties is explained by the fact that 
the inversion errors are determined by the error in the bulk of the frequencies (which are MDI 
intermediate-/ frequencies), and not just the errors in the low-degree modes. 

The MDI and GONG sets have good coverage of intermediate degree modes. The MDI set 
has p-modes from £ = up to a degree of £ = 194 while the GONG set has modes from £ = 
up to £ = 150. However, both these sets are somewhat deficient in low degree modes. The 
LOWL1 + BiSON combination, on the other hand, has a better coverage of low degree modes, but 
has modes from £ = only up to £ = 99. The GOLF data sets only contain low-degree modes 
(/ = 0, 1, 2) and hence have to be combined with other data before they can be used to determine 
the solar sound speed profile. 
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Fig. 14. — Five precise helioseismological measurements versus BP98. The figure compares the 
fractional difference between the sound speeds calculated for the 1998 Standard solar model (BP98) 
and the sound speeds in five helioseismological experiments. The references to the helioseismological 
data are given in the text. (The GOLF2 data were not available when this comparison was originally 
made.) The rms fractional difference between the calculated and the measured sound speeds is 
0.13% for all solar radii between between 0.05-Rq and O.95i?0. 



Figure 14 shows the somewhat less precise agreement that was obtained between BP98 (Bah- 
call, Basu, & Pinsonneault 1998) and the helioseismological data. The BP98 model and the observed 
sound speeds agree in the solar interior to about the same accuracy as for BP2000 and the observed 
speeds. However, the BP98 model sound speeds are about 0.1% smaller than the observed speeds 
in the broad intermediate region between 0.3R/Rq and 0.7R/Rq (cf. Figure 13 and Figure 14 ). 
Averaged over the entire region over which good measurements are available, 0.05-R© to O.95i?0, 
the rms fractional difference between the BP98 model and the LOWL1 + BiSON sound speeds is 
13 x 10~ 4 , which sound be compared with a rms difference of 10.4 x 10" 4 for the BP2000 model 
(see § 6.1) 5 . 



5 The Standard model described in this paper differs from the BP98 model in two respects, both discussed in § 2.1, 
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In BP98, we speculated that the broad feature of disagreement at the 0.1% level in Figure 14 
might be due to a combination of small errors in the adopted radiative opacities or in the equations 
of state. We investigated the implications for solar neutrino fluxes of the possibility that the small, 
broad discrepancy was due to an opacity error and concluded that if this were the case then the 
corrected solar model would predict 7 Be and 8 B neutrino fluxes that are about 5% larger than the 
fluxes predicted by BP98. 

Indeed, the origin of the broad discrepancy is due to an error in interpolating the radiative 
opacity near the edges of the opacity tables, as explained in § 2.1. The effect on the neutrino fluxes 
is somewhat smaller than we had estimated, an increase (relative to BP98) of 2% for 7 Be neutrinos 
and 3% for 8 B (cf. the results for GN93 in Table 7 and Table 1 of BP98). 

6.3. Sounds speeds: rotation 

Figure 15 compares the sound speeds of the Rotation model with the six precise observed sets of 
sound speeds. In the region between 0.3-R© and 0.6R®, the agreement between the Rotation model 
and the helioseismological data is slightly less good than with the Standard solar model (compare 
Figure 13 and Figure 15). Quantitatively, over the entire region between 0.05-R© to 0.95-Rq, the 
rms fractional difference between the Rotation model and the LOWL1 + BiSON sound speeds is 
12 x 1CT 4 , which is to be compared with 10 x 10~ 4 for BP2000. In the solar core (< O.25i? ), the 
two models, Rotation and BP2000, have almost identical rms fractional differences with respect to 
the LOWL1 + BiSON sound speeds, 0.073% and 0.0064%, respectively. 

that are significant for helioseismology: the correction of the opacity interpolation error and the slightly different 
heavy element mixture adopted for the BP2000 model. For helioseismology, the only significant different between 
GN93 and BP98 is the correction of the interpolation error; the heavy element abundances used for both models are 
the same (with Z/X — 0.0245). When a figure like Figure 13 was constructed using GN93 instead of BP2000, the 
general form of the small differences was practically the same. In fact, the rms difference averaged over the whole 
sun between GN93 and LOWL1 + BiSON is 0.00086, which is slightly better than the value of 0.00104 for BP2000 
versus LOWL1 + BiSON. One therefore obtains the correct impression by comparing Figure 13 and Figure 14. 
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Fig. 15. — Six precise helioseismological measurements versus the Rotation model. The figure com- 
pares the fractional difference between the sound speeds calculated for the Rotational solar model 
and the sound speeds in six helioseismological experiments. The model was developed by Pinson- 
neault and collaborators to explain the depletion of lithium. The references to the helioseismological 
data are given in the text. 



Figure 16 compares both the Rotation and the Standard model with just the LOWL1 + BiSON 
sound speeds, This figure shows that the Rotation model gives marginally better agreement with 
the measured sound speeds right at the base of the convective zone, comparable agreement in the 
deep interior, r < O.25i?0, and slightly less good agreement in the intermediate region between 
O.3i? and 0.6.R©. 

The neutrino fluxes calculated with the Rotation model lie well within the estimated errors in 
the Standard model fluxes, as can be seen easily by comparing the fluxes and the errors given in 
Table 7 and Table 10. The 7 Be flux for the Rotation model is 3% less than the Standard model 7 Be 
and the 8 B flux is 5% lower than the corresponding Standard model value. The Rotation model 
predicts a capture rate by 37 C1 that is 0.3 SNU less than the Standard model rate and a 71 Ga 
capture rate that is 2 SNU less than the Standard model rate. 
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Fig. 16. — BP2000 versus the Rotation model. The figure compares the fractional difference 
between the sound speeds determined from the LOWL1 + BiSON data with the sound speeds 
calculated for the BP2000 solar model and the Rotational solar model. The BP2000 model agrees 
slightly better with the measured sound speeds in the intermediate region between 0.3-R© and 
0.6^0. 

7. Helium abundance and depth of the convection zone 

Table 14 gives the calculated present-epoch values for the depth of the convective zone, R(CZ), 
and the surface helium abundance, Y s , for all 11 of the solar models considered in this paper. 

The observed values determined using measurements of p- modes are R(CZ) = (0.713 ± 
O.OOl)i?0 (Basu &: Anita 1995; for earlier work see Christensen-Dalsgaard, Gough, & Thompson 
1991) and Y s = 0.249 ± 0.003 (Basu & Anita 1997; see also Richard et al. 1996). 

The quoted errors are systematic; the statistical uncertainties are much smaller. For example, 
the cited uncertainty for Y s is designed to span the two values obtained when using two different 
equations of state (cf. Basu &; Anita 1997). There is no rigorous way of establishing a confi- 
dence level based upon agreement within, for example, one or two times the estimated systematic 
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Table 14: Helium abundance and the depth of the convective zone. The calculated present-day 
helium abundance on the surface and the convective zone depth are given for 11 solar models 
discussed in this paper. The quoted errors in the measured values of Y and R(CZ) represent best 
estimates of the systematic uncertainties but cannot be interpreted rigorously in terms of la or 3a 
errors. 



Standard 


Y s 


R(CZ) 


O lidbLKAcLL U. 




71 a 


NACRE 


0.244 


0.713 


ASOO 


0.239 


0.714 


GN93 


0.245 


0.712 


Pre-MS 


0.246 


0.713 


Rotation 


0.248 


0.714 


Radius _ 78 


0.245 


0.712 


Radium^sos 


0.245 


0.712 


No diffusion 


0.266 


0.726 


Old physics 


0.248 


0.712 


<s 34 = 


0.242 


0.715 


Mixed 


0.254 


0.732 


Measured 


0.249 ± 0.003 


0.713 ±0.001 



uncertainty. We point out that the recent reassessment of Z/X by Grevesse and Sauval (1998) 
resulted in a value of Z/X = 0.0230 which differs from the helioseismologically recommended value 
of Z/X = 0.0245 ± 0.0008 by twice the quoted uncertainty of the helioseismological determination. 

As a rule-of-thumb, we shall regard agreement within two times the quoted systematic un- 
certainty as 'satisfactory' and agreement within one times the quoted systematic uncertainty as 
'excellent'. 

Nine of the eleven solar models considered in this paper give excellent or satisfactory agreement 
with the observed depth of the convective zone. The only exceptions to the good agreement with 
the measured convective zone depth are the No Diffusion and the Mixed models, which are both 
strongly disfavored by the helioseismological measurements. 

The No Diffusion, 534 = 0, and ASOO models are the only ones that are not within twice 
the quoted uncertainty in the measured surface helium abundance. The No Diffusion model yields 
a surface helium abundance that is more than five times the quoted uncertainty away from the 
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helioseismological measurement, which very strongly disfavors the No Diffusion model. The surface 
helium abundance calculated for the ASOO model is also rather far from the helioseismological 
value, more than three times the quoted uncertainty in the helioseismological determination. This 
discrepancy should be examined more fully in future years as the abundance determinations are 
refined and the helioseismological determination of the helium abundance is repeated and the input 
data to the helioseismological analysis is varied over a wide range of allowed possibilities. 

8. Discussion and summary 

This paper provides new information about four topics: 1) the characteristics of the Standard 
solar model at the current epoch (see § 8.1 below); 2) time dependences of important characteristics 
of the Standard solar model (§ 8.2); 3) neutrino fluxes and related quantities for standard and 
variant solar models (§ 8.3); and 4) measured versus calculated solar sound speeds (§ 8.4). Extensive 
numerical data that are useful for applications are available at http://www.sns.ias.edu/~jnb . 

Just for fun, we provide our favorite list, our "top three," among the many disparate results 
presented in this paper. Our top three results are listed in § 8.5. 

8.1. Standard solar model: current epoch 

We present detailed numerical tabulations of the computed characteristics of our Standard 
solar model, which is defined and discussed in § 2. These tables include, as a function of the 
solar radius, the enclosed mass fraction, the temperature, mass density, electron number density, 
the pressure, and the luminosity fraction created in a given spherical shell, as well as the mass 
fractions of 1 H, 3 He, 4 He, 7 Be, 12 C, 14 N, and 16 0. Over the years, previous numerical versions of 
our Standard model have been used for a variety of purposes that range from comparisons with other 
stellar evolution codes, estimating the importance in the sun of newly considered physical effects, 
searching for possible instabilities in the sun, comparison with helioseismological measurements, 
and the calculation of processes (especially the MSW effect) that influence the propagation of solar 
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neutrinos. 

In the past, we have published in hard copy form increasingly more detailed and precise 
numerical tables of the characteristics of the solar interior. The capabilities of current calculations 
and the requirements of some of the most interesting applications have made complete hard copy 
publication no longer appropriate. We have therefore limited ourselves in § 2 to describing briefly 
the ingredients we use in calculating the current Standard model. We present the numerical results 
in exportable data files that are available at http://www.sns.ias.edu/~jnb . 

8.2. Standard solar model: time dependences 

For the first time in this series of papers, we have focused, especially in § 3, on details of the 
time dependence of important characteristics of the Standard solar model. 

The total luminosity in the Standard model increases by 48% from the zero-age main sequence 
stage to the present epoch. Over the same period, the effective temperature varies by only ±1.3%. 
These predictions constitute constraints on models for the evolution of the earth. 

The predicted time evolution of the solar luminosity is robust. Figure 2 shows that all solar 
models, even those models with deficient physics that are strongly disfavored by helioseismological 
measurements, predict essentially the same luminosity evolution. The average rms deviation of 
the deviant models, the Mixed, No diffusion, and S34 = models, from the standard solar model 
luminosity is only 1% over the history of the sun from 1 Gyr to the current epoch (see § 3.1 for 
more details). 

Table 4 presents the calculated large and small separations of the p-mode frequencies as a 
function of age for the standard solar model. 

We have presented in § 3 the time evolution of some of the principal physical quantities char- 
acterizing the solar core (the central temperature, density, pressure, and hydrogen mass fraction, 
as well as the fractions of the solar luminosity generated by different nuclear reactions). We also 
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present the evolution of important quantities at the base of the convective zone (radiative opacity, 
temperature, density, and pressure). We hope that these data and the scaling relations we have 
inferred will be sufficient to permit a future physical understanding of the time dependences using 
analytic and semi- analytic arguments. 

We find some simple scaling relations. For example, the solar luminosity, L & (t), is approxi- 
mately related to the solar radius, RQ,(t), as L@(t) oc Rg(t) 2 ' 5 . The depth of the convective zone, 
R(CZ,t), scales as i?(CZ,i) oc Ro(t), and the central temperature, T c (t), shows a similar behav- 
ior, T c (t) oc Ro(t). Moreover, we find that the mass of the convective zone, M(CZ,f) satisfies 
M(CZ,i) oc i? Q (i)~ 2 . The effective temperature is approximately constant, varying by only ±0.7% 
from a solar age of 2 billion years to 8 billion years. 

These results make predictions that are potentially testable. In principle, the measurement 
of the luminosity (by astrometry) of a star with the same mass and chemical composition (Z/X) 
as the sun would allow the prediction of the star's effective temperature, the depth and mass 
of the convective zone, and the large and small separations of the p-mode frequencies (see for 
example Monteiro, Christensen-Dalsgaard, & Thompson 2000 and Christensen-Dalsgaard 1997). 
In practice, it is difficult to make measurements sufficiently accurately to make possible precise 
tests of stellar evolution theory. For an appraisal of both the potential and the difficulty of making 
such measurements, the reader is referred to the recent papers by Guenther and Demarque (2000) 
and Morel, Provost, Lebreton, Thevenin, & Berthomieu (2000) on the binary pair of approximately 
solar mass stars, a Centauri AB (1.1M Q and O.9M ). 

8.3. Neutrino fluxes and related quantities 

Figure 8 and Table 11 give the electron number density as a function of position in the sun 
for the Standard solar model. The distribution of the electron density is required to compute 
the probability for matter-induced oscillations between active neutrinos. Similarly, Figure 9 and 
Table 12 give the radial distribution of the number density of scatterers of sterile neutrinos, n ste riie> 
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in the Standard solar model. We have not previously published precise values for the electron 
number density, or of n s t er iie> i n the outer regions of the sun. The outer regions are relevant for large 
mixing angle neutrino oscillations with relatively low neutrino mass differences (Am 2 < 10 _8 eV 2 ). 

Table 7 presents the neutrino fluxes and the event rates in the chlorine, gallium, lithium, and 
electron-scattering neutrino experiments that are predicted by the Standard solar model. These 
predictions assume that nothing happens to solar neutrinos after they are produced. The table also 
gives estimates of the uncertainties in the fluxes and the event rates; § 5.1 contains a discussion of 
the physical origin of the uncertainties, as well as the software used to calculate the asymmetric 
error estimates. 

How do the predictions of solar neutrino event rates compare with experiment? Table 8 
compares the predictions of BP2000 with the results of the chlorine, GALLEX + GNO, SAGE, 
Kamiokande, and Super-Kamiokande solar neutrino experiments. This table assumes nothing hap- 
pens to the neutrinos after they are created in the sun. The standard predictions differ from the 
observed rates by many standard deviations. Because of an accidental cancellation, the predicted 
solar neutrino event rates for BP2000 and BP98 are almost identical (see § 5.1.1). 

Table 10 compares the neutrino fluxes and the experimental event rates for all nine of the 
solar models whose helioseismological properties were investigated in BPB2000, plus two additional 
standard-like models considered here which have somewhat different heavy-element to hydrogen 
ratios. The seven standard-like models (the first seven models in Table 10) all produce essentially 
the same neutrino predictions; the spreads in the predicted pp, 7 Be, and 8 B fluxes are ±0.7%. 
±3%. and ±6.5% , respectively. The calculated rates for the seven standard-like solar models have 
a range of ±0.45 SNU for the chlorine experiment and ±2 SNU for the gallium experiments. 

The estimated total errors from external sources (see Table 7), such as nuclear cross section 
measurements and heavy element abundances, are about a factor of three larger than the uncer- 
tainties resulting from the solar model calculations. 

We have investigated one possible source of systematic errors, the relative weights assigned 
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to different determinations of nuclear fusion cross sections. We calculated the neutrino fluxes and 
predicted event rates using the NACRE (Angulo et al. 1999) fusion cross sections rather than the 
Adelberger et al. (1998) cross sections. The NACRE parameters lead to slightly higher predicted 
event rates in solar neutrino experiments. However, all changes in the neutrino fluxes and event 
rates between the NACRE-based predictions and the Standard predictions (based upon Adelberger 
et al. nuclear parameters) are much less than the la uncertainties quoted for the Standard model 
(cf. Table 7 and Table 9). 

The neutrino event rates predicted by all seven of the standard-like solar models considered 
here are inconsistent at the 5<r level (combined theoretical and experimental errors) with the results 
of the two gallium experiments, GALLEX and SAGE, assuming no new physics is occurring. The 
inconsistency with the chlorine experiment is similar but more complex to specify, since the largest 
part of the theoretical uncertainty in the calculated standard capture rate is due to the electron- 
type neutrinos from 8 B beta-decay. The fractional uncertainty in the 8 B flux depends upon the 
magnitude of the flux created in the solar interior. Moreover, the amount by which this 8 B flux is 
reduced depends upon the adopted particle physics scenario. 

A similar level of inconsistency persists even for the ad hoc deficient models, such as the S34 = 
and Mixed models, that were specially concocted to minimize the discrepancy with the neutrino 
measurements. For example, the calculated rates for the Mixed model are 6.15+Qg 5 SNU for the 
chlorine experiment and 115^' 8 SNU for gallium experiments. The Mixed model is 4.2a below the 
measured chlorine rate and 6.3c below the measured gallium rate. In addition, the deficient models 
are strongly disfavored by the helioseismology measurements. 

Figure 10 shows the calculated time dependence of the pp, 7 Be, 8 B, and 13 N solar neutrino 
fluxes. At the current epoch, the pp flux is increasing at a rate of 4% per billion years, and the 
7 Be, 8 B, and 13 N fluxes at the rates of 45%, 90%, and 99% per billion years, respectively. Since 
the age of the sun is estimated to be uncertain by only 5 x 10~ 3 billion years (see the Appendix by 
Wasserburg in Bahcall & Pinsonneault 1995), the age of the sun does not represent a significant 
uncertainty for solar neutrino predictions. 
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8.4. Sound speeds 

Figure 11 shows the excellent agreement between the helioseismologically determined sound 
speeds and the speeds that are calculated for the Standard solar model. The scale on this figure 
was chosen so as to highlight the contrast between the excellent agreement found with the Standard 
model and the two orders of magnitude larger rms difference for a solar model that could reduce 
significantly the solar neutrino problems. One would expect characteristically a 9% rms difference 
between the observations and the predictions of solar models that significantly reduce the conflicts 
between solar model measurements and solar model predictions. Averaged over the entire sun, the 
rms fractional difference is only 0.10% between the Standard solar model sound speeds and the 
helioseismologically-determined sound speeds. The agreement is even better, 0.06%, in the interior 
region in which the luminosity and the neutrinos are produced. Table 13 presents numerical values 
for the sound speeds predicted by the Standard solar model at representative radial positions in 
the sun. 

All eight of the standard-like solar models considered in this paper give acceptable agreement 
with the measured depth of the convective zone and the surface helium abundance (see Table 14). 
Of the twelve Standard, standard- like, variant, and deviant solar models considered here, only the 
No Diffusion and Mixed models disagree strongly with the measured convective zone depth. 

There are small but robust discrepancies between the measured and the calculated solar sound 
speeds. Figure 13 shows the fractional differences between the Standard model sound speeds and 
the speeds measured in each of six different determinations (using p-mode data from LOWL1, 
BiSON, GONG, MDI, GOLF, and GOLF2). The vertical scale for Figure 13 is expanded 21 times 
compared to the vertical scale of Figure 11. As can be seen from Figure 13, the differences between 
observed and measured sound speeds are comparable over much of the sun to the differences between 
different measurements of the sound speed, but there is a clear discrepancy near the base of the 
convective zone that is independent of which observational data set is used. There may also be a 
less prominent discrepancy near 0.2R Q . 
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The agreement between the BP2000 sound speeds and the measured values is improved over 
what was found earlier with the BP98 model. This improvement may be seen by comparing 
Figure 13 and Figure 14. The improvement is striking in the region between 0.3R & and 0.7i? Q 
and is due to the correction of an error in the interpolation algorithm for the radiative opacity (cf. 
discussion in § 2.1). 

The Rotation model includes a prescription for element mixing that was designed to explain 
the depletion of lithium. The calculated differences between the Rotation and the Standard models 
represent a reasonable upper limit to the effects that rotation, sufficient to explain lithium depletion, 
might produce. The overall agreement between the sound speeds of the Rotation model and the 
helioseismologically determined sound speeds is slightly worse, 0.12% rather than 0.10%, than the 
agreement obtained with the Standard solar model. The differences between the solar neutrino 
fluxes predicted between the Rotation model and the Standard model is typically about 0.3a u , 
where a v represents the uncertainties given in Table 7 for the Standard model neutrino fluxes. 
We conclude that further improvements of the theoretical calculations motivated by refinements of 
p— mode oscillation measurements are unlikely to significantly affect the calculated solar neutrino 
fluxes. 

8.5. Top three results 

Here are our favorite three results in this paper. 

• The robust luminosity evolution of the sun (see Figure 2). 

• The excellent agreement of the Standard model sound speeds with the measured sound speeds 
on the scale relevant for solar neutrino discussions (see Figure 11). 

• The simple relations as a function of time between the solar radius, the solar luminosity, the 
depth of the convective core, and the mass of the convective core (see Eq. 1-Eq. 5). 

We are grateful to many colleagues in solar physics, nuclear physics, and particle physics for 
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